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ABSTRACT 

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical 
CoRoT light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have 
infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian co- 
rotation radius. The repeating photometric dips have FWHM generally less than one day, depths 
almost always less than 15%, and periods (3 < P < 11 days) consistent with dust near the Keplerian 
co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually 
persist for several weeks and, in two cases, were present in data collected on successive years. For 
several of these stars, we also measure the photospheric rotation period and find that the rotation 
and dip periods are the same, as predicted by standard “disk-locking” models. We attribute these 
flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to 
the stellar photosphere. In some cases, these dips are also present in simultaneous Spitzer IRAC light 
curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars 
with light curves exhibiting this behavior to other classes of YSO in NGC 2264. A number of physical 
mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several 
of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with 
the stellar magnetic field or dust entrained in funnel-flow accretion columns arising near the inner 
disk wall. 

Subject headings: open clusters and associations: individual (NGC 2264) — circumstellar matter — stars: 
pre-main sequence — stars: protostars — stars: variables: T Tauri 


stauffer@ipac.caltech.edu 

* Based on data from the Spitzer and CoRoT missions, as 
well as the Canada France Hawaii Telescope (CFHT) MegaCam 
CCD, and the European Southern Observatory Very Large Tele- 
scope, Paranal Chile, under program 088.C-0239. The CoRoT 
space mission was developed and is operated by the French space 
agency CNES, with particpiation of ESA’s RSSD and Science 
Programmes, Austria, Belgium, Brazil, Germany, and Spain. 
MegaCam is a joint project of CFHT and CEA/DAPNIA, which 
is operated by the National Research Council (NRC) of Canada, 
the Institute National des Sciences de l’Univers of the Centre 
National de la Recherche Scientifique of France, and the Univer- 
sity of Hawaii. 

1 Spitzer Science Center, California Institute of Technology, 
Pasadena, CA 91125, USA 

2 Astronomy Department, California Institute of Technology, 
Pasadena, CA 91125, USA 

3 Jet Propulsion Laboratory, California Institute of Technol- 
ogy, Pasadena, CA 91109, USA 

4 Department of Physics and Astronomy, 5151 State Univer- 
sity Drive, California State University at Los Angeles, Los An- 
geles, CA 90032 

5 Centro de Astrobiologfa, Dpto. de Astroffsica, INTA-CSIC, 
PO BOX 78, E-28691, ESAC Campus, Villanueva de la Canada, 
Madrid, Spain 

6 Departamento de Ffsica - ICEx - UFMG, Av. Antonio Car- 
los, 6627, 30270-901, Belo Horizonte, MG, Brazil 

7 Universite de Grenoble, Institut de Planetologie et 
d’Astrophysique de Grenoble (IPAG), F-38000 Grenoble, France; 
CNRS, IPAG, F-38000 Grenoble, France 

8 Department of Astronomy, University of Michigan, 500 
Church Street, Ann Arbor, MI 48105, USA 

9 INAF - Osservatorio Astronomico di Palermo, Piazza del 
Parlamento 1, 90134, Palermo, Italy 


10 Department of Physics and Astronomy, The University of 
Georgia, Athens, GA 30602-2451, USA 

11 Department of Astronomy, University of Massachusetts, 
Amherst, MA 01003, USA 

12 U.S. Naval Observatory, Flagstaff Station, 10391 West 
Naval Observatory Road, Flagstaff, AZ 86001, USA 

13 Department of Physics and Astronomy (MS-9164), Western 
Washington Univ., 516 High St., Bellingham, WA 98225, USA 

14 NASA/Goddard Space Flight Center, Code 665, Greenbelt, 
MD 20771, USA 

15 Astronomy Department, Wesleyan University, Middletown, 
CT 06459, USA 

16 Department of Physics, University of Oxford, Keble Road, 
Oxford, OX1 3RH, UK 

17 European Space Agency, 8-10 rue Mario Nikis, F-75738 
Paris Cedex 15, France 

18 NASA Sagan Fellow 

19 Departamento de Ffsica e Matematica - UFSJ - Rodovia 
MG 443, KM 7, 36420-000, Ouro Branco, MG, Brazil 

20 NASA Ames Research Center, Kepler Science Office, Moun- 
tain View, CA 94035 


2 


1. INTRODUCTION 

Understanding how young stellar objects (YSOs) as- 
semble themselves and how the accretion process works 
requires understanding the structure and evolution of the 
inner disk on size scales < 0.1 AU. Even with the advent 
of infrared (IR) and radio interferometry, imaging of even 
the nearest YSOs to detect structure on these size scales 
is impractical. It therefore becomes necessary to infer the 
spatial structure and composition on these size scales in- 
directly from, for example, spectral energy distributions 
(SEDs), high- resolution spectra, or variability studies. 

SEDs have provided some of the most important in- 
sights into the structure and evolution of the inner disks 
of YSOs. These SEDs have been used not only to argue 
for the presence of an inner disk wall, located roughly 
at the dust sublimation boundary, but also that the up- 
per edge of the wall is in most cases curved (McClure 
et al. 2013). Some YSOs have both inner and outer op- 
tically thick disks separated by an optically thin gap, 
with the gap likely the result of gravitational sweeping by 
an embedded giant planet (Quillen et al. 2004). Mid-IR 
spectrophotometric variability seen in a number of YSOs 
suggest that the height of their inner disk wall varies by 
up to 20% on timescales of days or weeks (Espaillat et 
al. 2011). 

Optical or near-IR photometric variability studies offer 
another path to understanding inner disk structure and 
to the physics of the interaction between the inner disk 
and the parent star. Such studies date back more than 
half a century (Joy 1945), where variability was used 
as one of the defining characteristics of the new class 
of young, forming stars. Herbst et al. (1994) summa- 
rized the available information on T Tauri star variabil- 
ity, identifying several different classes dominated vari- 
ously by rotational modulation of cold or hot spots, vari- 
able accretion or variable extinction. Variable extinction 
is particularly attractive as a means to study inner disk 
structure because in principle it can yield information 
on very small size scales, if the cadence of the time series 
photometry is high enough. Bouvier et al. (1999) showed 
that AA Tau has broad, periodic flux dips whose ampli- 
tude and shape vary from epoch to epoch, which they 
ascribe to variable extinction resulting from a warp in 
the star’s inner disk and to large variation in the amount 
of material at the inner disk rim on ~week timescales. 
Photometric surveys of entire star forming regions (e.g., 
Carpenter et al. 2001; Flaherty et al. 2012; Findeisen et 
al. 2013) with modern 2-D arrays have also allowed other 
types of variable extinction events - both periodic and 
aperiodic - to be identified. 

We have recently obtained the highest cadence, high- 
est signal-to-noise (S/N) time domain survey of a large 
population of YSOs in a star-forming region during two 
campaigns in 2008 and 2011. We have used those data 
for NGC 2264 for an initial study of the AA Tau- type 
variables in the cluster (Alencar et al. 2010); to provide 
a statistical basis for a new YSO light-curve morphology 
classification scheme (Cody et al. 2014); and to iden- 
tify a set of YSOs whose light curves are dominated by 
short duration accretion bursts (Stauffer et al. 2014). In 
this paper, we use this dataset to identify a small set 
of YSOs whose light curves exhibit short duration, of- 
ten near-Gaussian shaped, periodic flux dips which must 


arise from dust structures near the inner disk wall. In 
a companion paper (McGinnis et al. 2015), we provide 
a quantitative analysis of all of the other YSOs in our 
NGC 2264 sample with flux dips due to variable extinc- 
tion events. 

2. OBSERVATIONAL DATA 

A detailed description of the photometric data we col- 
lected in our 2008 and 2011 observing campaigns for the 
NGC 2264 star-forming region is provided in Cody et 
al. (2014)p] The subset of the data sources that we use 
here include several sets of time-series, photometric mon- 
itoring data (the first four items below) plus additional 
single-epoch ancillary photometry (the last two items): 

* A 2008 campaign using the CoRoT (Convection, Rota- 
tion, and planetary Transits; Baglin et al. 2009) satellite, 
designated as SRaOl within the CoRoT mission, with 
photometry being collected from roughly MJD 54534 to 
54556, a period of 23 days. The observing cadence was 
generally every 512 seconds; however, bright targets and 
targets of known special interest could be done at 32 
second cadence (and the spacecraft could autonamously 
switch to 32 second cadence if transit-like events were 
detected). The passband is set just by the sensitivity 
curve for the CoRoT CCD, corresponding to a broad V 
+ R filter. Only stars present in an input catalog have 
CoRoT light curves. For each cataloged star, a specific 
pixel mask was designated to determine which set of pix- 
els was queried to derive the photometry for that star. 
For the light curves discussed in this paper and illus- 
trated in the figures, typical RMS uncertainties in the 
CoRoT relative photometry range from about 2 mmag 
at the bright end (e.g. Mon-250, Figure to about 6 
mmag at the faint end (Mon-6975, Figurejlj). 

* A 2011 campaign designated as SRa05 within the 
CoRoT mission, with data collected from roughly MJD 
55897 to 55936. The noise level in these light curves is 
both higher and more spatially variable due to increased 
CCD radiation damage. In some cases, the same star was 
observed with different pixel masks in the two epochs, 
which can also change the derived flux. The observing 
cadence and basic pipeline reduction steps were, however, 
the same for the two campaigns. 

* Spitzer (Werner et al. 2004) warm mission Infrared Ar- 
ray Camera (IRAC; Fazio et al. 2004) 3.6 and 4.5 /.tm 
data obtained simultaneous with the CoRoT SRa05 cam- 
paign, roughly from MJD 55900 to 55930. The cadence 
was approximately every 2 hours. Stars near the lower 
or upper edge of the mapped region have data in only 
one of the two channels. All of the Spitzer light curves 
discussed in this paper are for stars that are fairly bright 
at IRAC wavelengths, with expected RMS uncertainties 
of order 1% for their relative photometric accuracies. 

* /-band imaging using the US Naval Observatory 
(USNO) 40” telescope obtained on clear or partially clear 
nights during two campaigns (winter 2010 and winter 
2011), roughly MJD 55505 to 55560 and MJD 55890 to 
55990. Typically, 100-200 epochs of data were collected 

22 The CoRoT and Spitzer light curves for all probable 
NGC 2264 members, as well as our broad band photometry for 
these stars, are available at http : //irsa. ipac . caltech.edu/data/ 
SPITZER/CSI2264 
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in 2010, and of order 900 epochs in 2011/2012. For each 
campaign, imaging was obtained for an approximately 
23’x23’ region, centered on the cluster center. Aperture 
photometry was performed for all stars in all images. 
The zero-points for each CCD frame were established by 
identification of a set of non-variable stars in the region 
surveyed, and forcing the photometry for these stars to 
be constant across all images. Each light curve was nor- 
malized to a mean relative magnitude of zero. For the 
stars in Table [I] the USNO relative photometry should 
have RMS uncertainties of less than 1% at the bright 
end, rising to of order 2% for the faintest of our stars. 

* Single epoch ugri imaging of the entire cluster, with 
standard Sloan Digital Sky Survey (SDSS) filters, ob- 
tained using Megacam on the Canada-France-Hawaii 
Telescope (CFHT) during two observing runs in 2010 
and 2012. These data are presented in detail in Venuti 
et al. (2014). 

* Spitzer IRAC 3.6, 4.5, 5.8 and 8.0 /im and MIPS 24 
(im single epoch imaging photometry of the entire cluster 
obtained early in the cryogenic mission, and reported in 
Sung et al. (2009) and Teixeira et al. (2012). 

3. OVERVIEW OF YSOS WHOSE LIGHT CURVES ARE 
DOMINATED BY VARIABLE EXTINCTION EVENTS 

If a clump of circumstellar dust passes through our 
line of sight to a YSO, the light from the star will dim 
by an amount roughly proportional to the column den- 
sity of dust, modified according to the size distribution 
of the grains (if the clump is optically thin), or by an 
amount proportional to the fractional area of the star 
that is covered (if the clump is optically thick but small 
in angular size relative to the star). If the dust clump is 
long lived and relatively fixed in location, the flux dips 
it produces should be periodic at the Keplerian rotation 
period for that clump. For some stars, a warped inner 
disk has been invoked (e.g., AA Tau Bouvier et al. 
1999; Flaherty & Muzerolle 2010) to explain deep flux 
dips that extend over half or more of the star’s rotation 
period. Dust clumps produced by instability processes 
in the disk, and levitated above the disk plane into our 
line of sight, could produce flux dips that do not repeat 
(Turner et al. 2010). Other longer timescale variable ex- 
tinction events could be associated with Keplerian rota- 
tion in the outer disk bringing dust into our line of sight 
(Bouvier et al. 2013); with orbital motion of a binary 
system causing our line of sight to intersect a circumbi- 
nary disk (Plavchan et al. 2008a); or with dust entrained 
in a spatially or temporally variable disk wind (Bans & 
Konigl 2012). 

In our 40+ day monitoring of the star-forming region 
NGC2264 with CoRoT and Spitzer , we have identified 
three different morphological light curve signatures which 
appear to be associated with variable extinction events. 
Figure [T] shows light curves for NGC 2264 Classical T 
Tauri stars (CTTS) which illustrate each of these differ- 
ent modes. These three modes are: 

• AA Tau analogs (examples illustrated in Figure [l£i 
and b.) Light curves of this type have relatively 
flat continua, punctuated by deep, broad flux dips 
occurring at the Keplerian orbital period for mate- 
rial near the inner disk rim. The generally accepted 
physical explanation is that there is a warp in the 


inner disk, and that a portion of this warp occults 
our line of sight (Bouvier et al. 1999); 

• Aperiodic extinctors (examples illustrated in Fig- 
ure [lj: and d.) Light curves of this type have 
relatively flat continua and deep, narrow and/or 
broad flux dips but with no obvious periodicity. A 
possible physical mechanism is dust levitated rela- 
tively high above the disk mid-plane by magneto- 
resonance instabilities (MRI) or other instabilities, 
which then settles back towards the mid-plane, ac- 
cretes onto the star, or escapes or sublimates on 
timescales of hours to days. See McGinnis et al. 
(2015) for a discussion of these stars; 

• Short duration, periodic dippers (newly identified 
group - examples illustrated in Figure {Tp and f). 
Light curves of this type have relatively shallow, 
relatively narrow flux dips with periods generally 
in the range 3-10 days, and with dip shapes of- 
ten approximately Gaussian. Sometimes these dips 
are superposed on a relatively flat continuum, but 
sometimes the dips are superposed on sine-wave 
light variations which we attribute to other phys- 
ical mechanisms. A physical mechanism for these 
flux dips must also involve dust structures in or 
near the inner disk rim that are brought into our 
line of sight by orbital motion. We discuss options 
for the origin of these structures in §7. 

There are other types of variables - primarily eclipsing 
binaries (or systems with transiting exoplanets) - which 
also have narrow flux dips in their light curves, but their 
light curves are distinctively different from those illus- 
trated throughout this paper. Figure [2] shows two eclips- 
ing binary (EB) light curves from the 2008 NGC 2264 
campaign. Figure is the light curve for Mon-256, 
a dM-dM binary member of NGC 2264 (Gillen et al. 
2014) that is a borderline CTTS (broad and structured 
Ha profile; modest IR excess); this is the only EB in 
NGC 2264 with an IR excess presumably due to circum- 
stellar or circumbinary dust. For well-detached systems, 
eclipses by stellar-sized bodies yield flux dips that are 
of much shorter duration (with a Gaussian fit width of 
a = 0.045 days for the eclipses of Mon-256) than those 
which we study in this paper. The eclipse depths for 
Mon-256 appear to vary considerably with time primar- 
ily due to the varying continuum level - once that is 
removed, the depths are quite stable, though slightly dif- 
ferent at the primary and secondary eclipse since each 
star has slightly different radius and effective temper- 
ature. Figure [2b shows the light curve of a field star 
projected towards NGC 2264; its light curve is that for a 
near-contact binary (P ~ 11 day) where the continuum 
variations are primarily due to the ellipsoidal shapes of 
the stars superposed on which are narrow flux dips from 
the primary and secondary eclipses. The Gaussian width 
for the eclipse dip in this case (a = 0.21 days) is within 
the range of the widths for the short-duration flux dips 
we discuss in this paper. However, the regularity imposed 
in an EB system yields a very stable light curve which is 
easily distinguishable from light curves associated with 
variable circumstellar extinction. 
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Fig. 1. Light curves for six CTTS members of NGC 2264 whose CoRoT data are dominated by flux dips due to variable extinction. The 
top two are AA Tau analogs, with broad flux dips most naturally explained as due to a warped inner disk intersecting our line of sight; the 
middle two light curves have aperiodic, narrow flux dips, possibly associated with MRI (magneto- resonance instability) events levitating 
dust in the inner disk; the bottom pair of plots exemplify the class of light curves with periodic narrow dips which are the subject of this 
paper (red dots provided to guide the eye to the periodic flux dips). Time zero in all panels except (e) corresponds to MJD = 55896.0 
days; for panel (e) time zero corresponds to MJD= 54533.0 days. Note, for Mon-21, the light curve shown has been processed to remove 
a smoothly varying “continuum” which we ascribe to photospheric spots - see Section 3.1.1. Compared to the AA Tau analogs (panels a 
and b), the objects of this study (panels e and f) show, on average, narrower, shallower and more symmetrical flux dips (see Figure |3|). 


CoRoT mag 


5 




Fig. 2. — CoRoT light curves for a) Mon-256 ( CoRoT SRaOl 223992193), an M dwarf eclipsing binary that is a member of NGC 2264; 
and b) Mon-6385 ( CoRoT SRaOl 223958210), a much higher mass field star that is a near-contact binary. The flux dips produced by 
eclipses in such systems are easily distinguished from those due to the variable extinction events we study here. 
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4. IDENTIFICATION AND INITIAL CHARACTERIZATION 
OF YSOS WITH SHORT-DURATION, PERIODIC FLUX 
DIPS IN THEIR LIGHT CURVES 

First results from our 2011-2012 intensive monitoring 
of NGC 2264 with CoRoT and Spitzer have been pub- 
lished in Cody et al. (2014) and Stauffer et al. (2014). 
While working on those papers, we noticed a small set 
of stars that shared some similarities but did not seem 
to fit into the primary groupings identified by Cody et 
al. Often, the “narrow” flux dips in these light curves 
are superposed on other, larger amplitude, types of vari- 
ability - primarily semi-sinusoidal light curves normally 
attributed to photospheric spots - making identification 
of this new light curve signature more difficult. 

These narrow flux dips presumably arise because some 
physical structure transits our line of sight to the stellar 
photosphere. As mentioned in the previous section, there 
is already a light curve class (AA Tau analogs) with peri- 
odic, variable extinction events and for which a standard 
physical model has been developed. It is possible that 
this new light curve class can be fit into the same physi- 
cal model; however, we believe there are a number of rea- 
sons to consider alternate physical models which we will 
discuss in §7. To emphasize the primary morphological 
differences between the “classical” AA Tau analogs and 
our proposed new class, Figure [3] compares the average 
widths and depths of the flux dips for the new group of 
stars to those of the A A Tau analogs in NGC 2264 identi- 
fied in the companion paper by McGinnis et alj^] While 
the distribution of points in Figure [3] could be regarded 
as simply one population with a range of properties, it 
could also be viewed as two populations whose properties 
only slightly overlap. The goal of this paper is to better 
characterize the light curves with these narrow flux dips, 
thereby determining if they should be considered as a 
new class of event or if they should simply be regarded 
as an extension of the AA Tau class. Throughout the 
paper, we use the terms “narrow” and “short-duration” 
interchangeably to describe the flux-dips of this set of 
stars, and “broad” (or long-duration) to describe the flux 
dips for the classic AA Tau light curves, with the divid- 
ing line being approximately FWHM/Period = 0.25 (see 
Figure [3]). 

In Tableland Table [2j we provide basic photometric 
and spectroscopic data for the nine stars which we believe 
fall into this “narrow-dip” group, compared to a total of 
159 NGC 2264 CTTS with CoRoT light curves (McGin- 
nis et al. 2015). All but one of these stars are previously 
identified members of NGC 2264 that had been classified 
as CTTS based on either Ha emission or SEDs. Spectral 
types are available for all of these objects, and all are 
low mass stars with spectral types of K5 or later. In the 
remainder of this section, we provide a brief description 
of each star and its light curve properties. In Table 0J 
we provide the fitted times, widths and depths of all the 
flux dips for these stars where those quantities could be 

23 Primarily because our definition of the AA Tau class requires 
that the optical light curves satisfy two quantitative tests (Cody 
et al. 2014), our (Stauffer and Cody) list of NGC 2264 AA Tau 
analogs is slightly smaller than that of McGinnis et al. (2015). 
Specifically, we do not include Mon-358, Mon-654, Mon-1054, and 
Mon-1167; we also do not include Mon-56 and Mon-1131 because 
we instead place them in our new category of stars with narrow, 
periodic flux dips. 


determined accurately, and we give a label for each of 
these dips - the dip ID column of the table - primarily as 
a means for us to refer to specific dips in the text. 

4.1. Mon- 21 

Mon-21, also known as V602 Mon, is a relatively poorly 
studied member of NGC 2264 of spectral type K5 (Rebull 
et al. 2002), with weak Ha emission (Furesz et al. 2006) 
and a relatively modest IR excess. It has no discernible 
ultraviolet (UV) excess (Venuti et al. 2014). We obtained 
C 0 R 0 Flight curves for it in both 2008 and 2011, and an 
IRAC light curve for it in 2011. 

4.1.1. Mon-21 (2008 CoRoT campaign) 

The as-observed C 0 R 0 Flight curve for Mon-21 in 2008 
(Figure |4^,) appears relatively complex because, we be- 
lieve, it is a superposition of two distinctly different wave- 
forms arising from two different physical mechanisms. 
Specifically, we ascribe the broad, semi-sinusoidal vari- 
ation to rotational modulation due to the presence of a 
large, high latitude spot group on the star’s photosphere. 
The narrow flux dips that are also present are the compo- 
nent we wish to isolate and quantify. The two waveforms 
appear to have nearly identical periods, and to produce 
features of very similar amplitude. 

In order to measure accurately the properties of the 
flux dips, it is necessary to remove the comparatively 
smoothly varying, semi-sinusoidal pattern. We construct 
a model light curve to emulate the sinusoidal variation by 
taking the CoRoT data from day MJD 4537.4 to 4540.6, 
the portion of the light curve least affected by narrow 
flux dips, and use those data to produce a sine-wave- 
like template covering one full period (3.22 days). (The 
light curve interval from ~4550 to 4553 is also relatively 
clean, but we felt that the inflection in the curve pre- 
sumably due to a narrow dip was harder to eliminate 
unambiguously compared to the earlier time interval.) 
From examination of the light curve, we mark beginning 
points (flux and time) for each cycle, which are indi- 
cated as the lower green circles in Figure 0b. A linear 
correction term is added to the flux for each replicated 
template so that its ending flux matches the next cycle’s 
beginning flux. Finally, we slightly adjust the amplitude 
of each waveform (by multiplying all the measurements 
in that replicated waveform by a constant) so that the 
maximum brightness matches the observed maximum at 
that epoch (upper green circles in Figure®)). The model 
light curve is shown in red in Figure 0b; the result of sub- 
tracting this model from the observed light curve - and 
hence our representation of the light curve solely due to 
variable extinction events - is shown in Figure |4]c. We 
made no further adjustments to the model to better fit 
the observed data other than the deterministic steps de- 
scribed above. We note that the first and last of the 
narrow dips are in fact double-dips, with the second dip 
being shallower and offset by of order 0.8 days in both 
cases (see Figure®). 

We have ascribed the sinusoidal variations to non- 
axisymmetrically distributed star spots on the stellar 
photosphere - probably cold spots. Could it instead 
be ascribed to variable extinction due to a warped disk 
intersecting our line of sight? We offer two arguments 
in favor of our hypothesis and counter to the alterna- 
tive. First, we know what the light curves of YSOs 
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Fig. 3. — Comparison of the average full depth of the flux dips to the ratio of the dip FWHM and the dip period for CoRoT light curves 
of “classical” AA Tau stars in our data set and for the stars of Table The stars in the newly identified narrow-dip class have flux dips 
that are both shallow and of relatively short duration. The stars marked with green plus signs represent AA Tau type variables also having 
periodic narrow flux dips, while the red plus signs represent s tars whose Co RoT light curves have flux dips that are sometimes narrow and 
sometimes broad, as discussed in §5 and illustrated in Figure |17| and Figure pfs] 


TABLE 1 

Basic Information for YSOs with Short Duration Flux Dips 


Mon ID a 

2MASS ID 

CoRoT 10 

CoRoT 

SpT d 

Ha EW 

(A) 

d FR4.5 e 

Period f 

(days) 

L„ g 

(L©) 

R* g 

(R©) 

Rco h 

(Rq) 

Rd 1 

(R©) 

CSIMon-000021 

06405944+0959454 

223980412 

223980412 

K5 

-7.4 

1.6 

3.2 

0.96 

1.69 

9.7 

7.3/14.5 

CSIMon-000056 

06415315+0950474 

223994760 

223994760 

K5 

-1.8 

1.0 

5.83 

0.77 

1.51 

14.3 

6.5/13.0 

CSIMon-000378 

06405292+0944544 


616872605 

K5.5 

-8.5 

2.1 

11.1 

0.90 

1.71 

21.3 

7.0/14.0 

CSIMon-000566 

06405275+0943004 

400007955 

400007955 

M3. 5 

-19.4 

0.9 

4.5 

0.16 

1.23 

7.3 

2. 9/5.9 

CSIMon-001076 

06403819+0929524 

400007702 

605424384 

Ml 

-2.6 

0.6 

5.85 

0.26 

1.25 

10.4 

3. 8/7.5 

CSIMon-001131 

06393441+0954512 

223957455 

223957455 

Ml. 5: 

-31.1 

1.4 

5.2 

0.55 

2.02 

8.9 

5.5/11.0 

CSIMon-001165 

06404192+0943551 



M3: 

-7.9 

1.3 

5.5 

0.13 

1.15 

8.5 

2. 7/5.3 

CSIMon-001580 

06410285+1007119 

223981250 


Ml: 

-17.2 

0.6 

7.1 

0.27 

1.28 

11.9 

3. 8/7.7 

CSIMon-006975 

06395318+1001137 


616803898 

M2. 5: 

-15. 

0.2 

2.8 

3.42 

3.44 

5.8 

9.3/18.6 

a The CSIMon IDs 

are our internal naming 

scheme for stars in the field of NGC 2264 — see 

Cody et al. 

(2014) and http://: 

irsa. ipac 

. caltech 

. edu/ 


data/SPITZER/CSI2264f orRAsandDecsandalternativelDs In the main text of this paper, we omit the “CSI” and the leading zeros in the object names 
for brevity. 

k CoRoT identification number in 2008 SRaOl campaign. 
c CoRoT identification number in 2011 SRa05 campaign. 

^ See Cody et al. (2014) and the Appendix for the sources of the spectral type and Ha equivalent width (EW) data. Here we use the convention 
that negative values are in emission, which is opposite to the convention in Cody et al. 

e FR4.5 is an estimate of the ratio of the flux from the disk to the flux from the stellar photosphere at 4.5 ^m, based on all the available broad-band 


photometry and the published spectral types. Because the SED models do not include veiling, these values are likely lower limits. 

Period determined in this paper, see Section 3. 
g Stellar luminosity and radius estimates from Venuti et al. (2014). 
k Co-rotation radius; see formula in Bouvier et al. (1999). 

1 Dust sublimation radius; see formula in Monnier 8z Millan-Gabet (2002). We adopt values of 1500 K for the dust sublimation temperature, and 
minimum and maximum values of their Qjj parameter of 1.0 and 4.0, based on their Figure 2. 


TABLE 2 

Photometric Data for the Stars of Table U 


Mon ID 

u 

9 

r 

i 

J 

H 

K a 

[3.6] 

[4.5] 

[5.8] 

[8.0] 

[24] 

CSIMon-000021 

18.174 

15.835 

14.810 

14.452 

12.781 

12.093 

11.685 

11.248 

10.936 

10.707 

10.052 

7.556 

CSIMon-000056 

18.785 

16.205 

15.056 

14.475 

12.959 

12.302 

12.003 

11.717 

11.390 

11.038 

10.290 

7.866 

CSIMon-000378 

18.277 

16.142 

14.943 

14.239 

12.699 

11.862 

11.267 

10.882 

10.540 

10.091 

9.371 

6.270 

CSIMon-000566 

20.988 

19.109 

17.747 

16.418 

14.318 

13.588 

13.279 

12.823 

12.576 

12.140 

11.132 


CSIMon-001076 

20.754 

18.336 

16.896 

16.012 

13.931 

13.194 

12.891 

12.485 

12.289 

12.025 

11.278 

8.566 

CSIMon-001131 

18.769 

17.180 

15.832 

14.856 

12.991 

12.215 

11.932 

11.492 

11.152 

10.893 

10.195 

6.998 

CSIMon-001165 

21.924 

19.528 

17.983 

16.427 

14.679 

13.845 

13.582 

13.043 

12.789 

12.380 

11.571 


CSIMon-001580 

20.489 

17.976 

16.553 

15.678 

13.978 

13.218 

12.874 

12.624 

12.455 

12.223 

11.244 

8.013 

CSIMon-006975 


20.027 

17.974 

16.626 

13.476 

12.576 

12.151 

11.666 

11.56 


11.10 



Note. — Broad-band photometry for the stars from Table [T| in AB magnitudes for ugri but in Vega magnitudes for longer wavelengths. The 
ugri data are from CFHT, as reported in Venuti et al. (2014); the JHK S data are from the on-line 2MASS all-sky point source catalog; the IRAC 
data are from Sung et al. (2009), or from our own analysis of archival IRAC imaging. Typical photon- noise and calibration uncertainties for these 
magnitudes are of order 0.02 mag; however, because all of these stars are photometric variables with amplitudes up to several tenths of a magnitude, 
these single-epoch data could differ from absolute, time-averaged values by 0.1 mag or more. 
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TABLE 3 

Flux Dip Gaussian Fit Parameters for Stars in Table |7] 


Mon ID 

Year 

Dip ID 

Center (days) 
MJD-50000. 

Width (days) 

Percent Depth 

CSIMon-000021 

2008 

A1 

4535.39 

0.17 

10% 

CSIMon-000021 

2008 

A2 

4536.15 

0.17 

7% 

CSIMon-000021 

2008 

B 

4541.99 

0.18 

7% 

CSIMon-000021 

2008 

C 

4545.07 

0.18 

5% 

CSIMon-000021 

2008 

D 

4548.47 

0.17 

12% 

CSIMon-000021 

2008 

El 

4554.73 

0.16 

10% 

CSIMon-000021 

2008 

E2 

4555.60 

0.15 

4% 

CSIMon-000021 

2011 

A 

5897.65 

0.22 

10% 

CSIMon-000021 

2011 

B1 

5900.35 

0.35 

2.5% 

CSIMon-000021 

2011 

B2 

5900.77 

0.29 

4.5% 

CSIMon-000021 

2011 

C 

5910.29 

0.26 

8% 

CSIMon-000021 

2011 

D1 

5913.10 

0.33 

3.5% 

CSIMon-000021 

2011 

D2 

5913.30 

0.11 

3.5% 

CSIMon-000021 

2011 

D3 

5913.64 

0.18 

6% 

CSIMon-000021 

2011 

E 

5916.31 

0.41 

4% 

CSIMon-000021 

2011 

FI 

5919.23 

0.38 

4.5% 

CSIMon-000021 

2011 

F2 

5919.73 

0.25 

6.5% 

CSIMon-000021 

2011 

Cl 

5922.51 

0.25 

6.0% 

CSIMon-000021 

2011 

G2 

5922.96 

0.20 

8.2% 

CSIMon-000021 

2011 

HI 

5925.04 

0.34 

1.5% 

CSIMon-000021 

2011 

H2 

5925.86 

0.29 

10% 

CSIMon-000021 

2011 

11 

5928.60 

0.35 

5.5% 

CSIMon-000021 

2011 

12 

5929.20 

0.27 

10% 

CSIMon-000021 

2011 

J1 

5931.79 

0.39 

4% 

CSIMon-000021 

2011 

J2 

5932.37 

0.23 

12% 

CSIMon-000056 

2011 

A 

5912.63 

0.32 

12.5% 

CSIMon-000056 

2011 

B 

5919.00 

0.25 

5.5% 

CSIMon-000056 

2011 

Cl 

5924.30 

0.25 

8.3% 

CSIMon-000056 

2011 

C2 

5925.08 

0.27 

7.7% 

CSIMon-001165 

2011 

A 

5911.87 

0.20 

25% 

CSIMon-001165 

2011 

B 

5922.83 

0.15 

14% 

CSIMon-006975 

2011 

A1 

5905.54 

0.10 

15% 

CSIMon-006975 

2011 

A2 

5905.82 

0.10 

12% 

CSIMon-006975 

2011 

B 

5908.36 

0.13 

14% 

CSIMon-006975 

2011 

C 

5910.94 

0.10 

7% 

CSIMon-006975 

2011 

D 

5919.59 

0.11 

11% 

CSIMon-001 131 

2011 

A 

5901.06 

0.30 

24% 

CSIMon-00 1131 

2011 

B 

5906.14 

0.30 

23% 

CSIMon-001131 

2011 

Cl 

5910.66 

0.27 

12% 

CSIMon-001 131 

2011 

C2 

5911.03 

0.18 

10% 

CSIMon-001131 

2011 

C3 

5911.56 

0.36 

11% 

CSIMon-001131 

2011 

D 

5926.75 

0.75 

13% 

CSIMon-001131 

2011 

El 

5931.66 

0.34 

6.5% 

CSIMon-001131 

2011 

E2 

5932.04 

0.14 

10% 

CSIMon-001131 

2011 

E3 

5932.31 

0.15 

12% 

CSIMon-001580 

2008 

A1 

4541.73 

0.18 

13% 

CSIMon-001580 

2008 

A2 

4542.15 

0.24 

9% 

CSIMon-001580 

2008 

B1 

4548.53 

0.15 

6% 

CSIMon-001580 

2008 

B2 

4548.7 

0.06 

5% 

CSIMon-001580 

2008 

B3 

4548.98 

0.17 

8% 

CSIMon-001580 

2008 

C 

4554.83 

0.20 

25% 

CSIMon-000378 

2011 

A 

5918.64 

0.55 

14% 

CSIMon-001076 

2008 

A 

4538.65 

0.17 

15.7% 

CSIMon-001076 

2008 

B 

4550.54 

0.16 

13% 

CSIMon-001076 

2008 

C 

4556.20 

0.15 

12.5% 

CSIMon-000566 

2011 

A 

5909.33 

0.32 

5% 
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MJD - 50000. (days) 


Fig. 4. — a) As observed CoRoT light curve for Mon-21; b) same as (a), except now overplotted in red by a model light curve designed 
to approximate the light curve shape attributable to spots on the star’s photosphere; c) residual light curve produced by subtracting the 
spotted star model from the observed light curve for Mon-21; d) superposition of the flux dips approximately centered at MJD 54536 (in 
blue, dips A1 and A2 in Table |3| and 54555 (dips El and E2, in red), highlighting their double-peak nature and their qualitative similarity; 
e and f) illustration that some of the narrow flux dips for Mon-21 are nearly Gaussian in shape (dips B and D in Table |3) respectively). 
The period corresponding to the spotted-star waveform - the red curve in panel (b) - is about 3.15 days; the period corresponding to the 
average dip interval is about 3.22 days (see text). 
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with typical warped disks look like because there are 
many of them in NGC 2264 (see the companion pa- 
per to this, McGinnis et al. 2015). Those light curves 
look nothing like the semi-sinusoidal pattern shown in 
red in Figure the warped disk light curves are 
much more irregularly shaped and much more time vari- 
able. Second, we also know that light curves dominated 
by non-axisymmetrically distributed, high-latitude spots 
are common amongst YSOs (Herbst 1986; Bouvier & 
Bertout 1989) and their light curves closely resemble the 
red curve in Figure |4|). In fact, we have identified one 
CTTS in NGC 2264 - Mon-103 - whose CoRoT light 
curve is nearly identical, except for minor scale factors, 
to the red curve in Figurepb. We provide a quantitative 
comparison of the 2008 Mon-21 light curve and that for 
Mon- 103 in the Appendix. 

Figures H}.’ and f compare a Gaussian fit to the observa- 
tions for the dips that are narrowest and best fit by such 
a model (dips B and D - see Tablep]) . The derived times, 
Gaussian widths and dip depths (as a percentage of the 
“continuum” flux) are provided in Table [3j Most of the 
dips are well fit, and most have widths m the narrow 
range from 0.15 to 0.18 days. 

The semi-sinusoidal pattern and the narrow flux dips 
both have apparent periodicities of about 3.2 days. How- 
ever, are the dips exactly periodic? And, are the periods 
for the semi-sinusoidal waveform and the dips precisely 
the same? The answer to the first question is no, as can 
be seen from examination of the data in Table 0 The 
centroids of each dip can be measured to an accuracy of 
a few hundredths of a day, yet dip B and C are separated 
by 3.08 days, while dips C and D are separated by 3.40 
days. Applying the period-finding algorithm of Plavchan 
et al. (2008b) to the light curve for just the flux dips (Fig- 
urefij:) yields a period of 3.224 days; a similar estimate of 
3.2213 days is obtained by subtracting the measured times 
for dip El and Al and assuming they are separated by 
six periods. The Lomb-Scargle algorithm yields a period 
of 3.15 days for the semi-sinusoidal pattern (represented 
either just by the red curve in Figure |4b or by the resid- 
ual remaining after subtracting the Gaussian fits from 
the original light curve) . Given the real variations in the 
timing of the flux dips and the comparative shortness of 
our CoRoT campaign, we cannot conclude the two peri- 
ods are significantly different. 

4.1.2. Mon-21 (2011 CoRoT campaign) 

By 2011, the CoRoT light curve for Mon-21 changed 
morphology, as illustrated in Figure [5k. The semi- 
sinusoidal pattern is either essentially absent or much 
reduced in amplitude. The flux dips appear somewhat 
broader and in many cases asymmetric. One possibility 
is that we are still seeing a blend of a spotted waveform 
and narrow flux dips from variable extinction, with the 
spot being responsible for the weaker feature causing the 
apparent asymmetry in each observed waveform; another 
possibility is that spots no longer contribute in any sig- 
nificant way and the observed waveforms are due to the 
blend of two narrow dips. In order to document the full 
light curve shape, we have chosen to simply model the 
observed features with Gaussian fits, usually with two 
components; the fits to a number of the dips are shown 
in Figure [5j On average, the primary dip at each epoch 
has a Gaussian sigma width of 0.26±0.02 days, about 


50% broader than the dips in 2008. The weaker sec- 
ondary dips always precede the deeper dip. See Table [3] 
for the fit parameters. The average period inferred from 
the dip centroids is 3.16±0.02 days, somewhat different 
from the dip period in 2008, but essentially equal to the 
period for the sinusoidal waveform we attribute to spots 
on the stellar surface. 

A ground-based USNO /-band light curve for Mon-21 
is available that overlaps in time with the CoRoT light 
curve. For a star of Mon-21’s I magnitude, the expected 
RMS uncertainty in the relative USNO photometry is 
about 0.02 mag. Because the CoRoT data have signifi- 
cantly lower noise and higher cadence, they provide the 
best measure of the shape of the flux dips. However, 
the USNO data extend beyond the CoRoT campaign 
by nearly two months, and therefore provide valuable 
information as to the persistence and amplitude vari- 
ability of the flux dips. Figure [BJi shows an expanded 
view of the region of best overlap between the CoRoT 
and USNO data; Figure^ shows the entire USNO light 
curve. These data illustrate that the same flux dip struc- 
ture is detected in the ground-based data, and that it per- 
sists for at least twenty days after the end of the CoRoT 
campaign. Because the ground-observing only provided 
data for about 5-6 hours per night, the window function 
for detecting the flux dips in the USNO data is not very 
good and it is not possible to deduce whether similar 
depth flux dips persisted beyond MJD 55960. 

We also obtained a Spitzer / IRAC [3.6] light curve for 
Mon-21 in 2011. As shown in Cody et al. (2014), for 
a star of Mon-21’s brighness, the expected RMS uncer- 
tainty in the relative IRAC photometry should be of or- 
der 1%. The most obvious feature of the IRAC light 
curve is a gradual flux decrease of 0.25 mag beginning 
around MJD 55910. On small timescales, there are pos- 
sible correspondences between the CoRoT flux dips and 
features in the IRAC data, but there are also places 
where there is no correspondence. Given the shallowness 
of the dips in the CoRoT data, and that the disk con- 
tributes equal or more light at [3.6], the lack of distinct 
dips in the IR is not unexpected. 

4.2. Mon-56 

In many respects, Mon-56 is a near twin to Mon-21. 
It also has a published spectral type of K5 (Dahm & Si- 
mon 2005), Ha emission with an inverse P-Cygni profile 
but quite small equivalent width (Furesz et al. 2006), es- 
sentially no UV excess and IR colors placing it near the 
center of the Class II population in the IRAC color-color 
diagram (Allen et al. 2004). Similar to Mon-21, we inter- 
pret its CoRo T light curve in 2011 as also being the su- 
perposition of two components with essentially identical 
periods - a smoothly varying, semi-sinusoidal component, 
possibly due to star-spots and a system of periodic, shal- 
low, short duration flux dips which are sometimes single 
and sometimes double. Figure [TJi shows the as-observed 
CoRo T light curve. 

Compared to Mon-21, however, the Mon-56 light curve 
is more difficult to deconvolve because the flux dips are 
generally stronger and broader and in phase with the 
flux minima of the sinusoidal waveform, thereby making 
it harder to construct a high quality model for the spot- 
ted light curve (or to be certain that this component is 
indeed due to spots as opposed to a warped-disk model). 
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MJD - 50000. (days) 



MJD — 50000. (days) 



MJD— 50000. (days) 




MJD -50000. (days) MJD-50000.(days) 


Fig. 5. — a) CoRoT light curve for Mon-21 in 2011 (in black). The blue crosses show the 2011 IRAC light curve for Mon-21, shifted in 
zero point to overlap the CoRoT data. The two light curves track each other fairly well at the start of the campaign and some of the later 
CoRoT flux dips have possible IRAC counterparts, but the 0.3 magnitude dimming at IRAC wavelengths beginning at MJD 55910 is not 
seen at optical wavelengths; b) flux dip C (see Table |3} which is reasonably well fit by a single Gaussian; c) flux dip H, where a shoulder to 
the flux dip requires the presence of a second narrow flux dip; d) flux dip J, also requiring a two component fit; and e) flux dip D, which 
appears to be well fit by a three Gaussian model. 
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Fig. 6. — a) A portion of the USNO I - band light curve for Mon-21, illustrating generally good agreement between the USNO and CoRoT 
data. The diamond symbols are plotted at the times of the deepest flux dips in the CoRoT data, b) The entire USNO /-band light curve. 
Additional flux dips appear to be present in the date range 55950 to 55960. 
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Fig. 7. — a) CoRoT light curve for Mon-56 in 2011; b) CoRoT light curve for Mon-56, with a model light curve designed to approximate 
the light curve shape attributable to spots on the star’s photosphere overplotted in red; c) Residual light curve for Mon-56, produced by 
subtracting the spotted star model from the observed light curve; d and e) CoRoT light curves for two of the prominent flux dips in Mon-56 
(dips A and B in Table |3|. The adopted perod for both the spotted-star waveform and the flux dips is about 5.8 days. 
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We have therefore resorted to a lower fidelity but sim- 
pler model for the smoothly varying component of the 
light curve - specifically, we adopt a simple sine- wave, 
whose phase, amplitude and period we set to yield the 
best fit (i.e. to yield a residual light curve with a nearly 
flat continuum level interspersed with narrow dips). Fig- 
ure I3> shows our final model for the smoothly varying 
component of the light curve superposed on the observed 
Mon-56 light curve, and Figure [7}; shows the result of 
subtracting this model from the observations, yielding 
our best visualization of the short-duration flux dips in 
this star. The full amplitude of the adopted sinusoid is 
0.06 mag (about 6%), and the period is 5.83 days (the 
same as determined from ground-based monitoring by 
Lamm et al. 2004). Because of the approximate nature 
of the deconvolution, we attempt to derive flux dip fits 
only for the three strongest features in Figure |7}:. Those 
parameters are provided in Table [3] and the fits are il- 
lustrated for the apparently single mix dips in Figure [7}l 
and e. 

We note that the Mon-56 CoRoT data are interpreted 
differently in McGinnis et al. (2015). McGinnis et al. at- 
tribute all of the variability seen in the system to variable 
extinction, as opposed to a superposition of spot and ex- 
tinction variability. Regardless of these differences, both 
models agree that the narrow, periodic flux dips are best 
attributed to extinction and we therefore believe Mon-56 
fits within the class of variability discussed in this paper. 

4 . 3 . Mon-1165 

Mon-1165 has spectral type M3 (see appendix). In an 
IRAC color-color diagram, it falls within the box defining 
the locus for Class II YSOs (Allen et al. 2004). However, 
as for the previous two YSOs, it has essentially no UV 
excess, and only a weak Ha emission line. 

We do not have a CoRoT light curve for Mon- 11 65, 
but the IRAC fight curve shows distinct, narrow, peri- 
odic, short-duration, paired flux dips (see Figurepm). An 
auto-correlation analysis of the light curve yields a pe- 
riod for the flux dips of 5.5 days. The lower amplitude 
dip always precedes the main dip by about 1.2 days (see 
Figure[8p). The flux dips are generally deeper at 3.6 pm 
than at 4.5 pm by, on average, about 50% but this is 
poorly determined due to the shallowness and narrow- 
ness of the dips and the relatively low cadence of the ob- 
servations. Only two of the dips are well-enough defined 
to accurately measure their shape and width - those fits 
are shown in Figure [8jc and d; the Gaussian er values for 
the two fits are 0.15 and 0.20 days. 

We have ground-based, optical, synoptic photometry 
which only overlap with a small portion of the Spitzer 
data due to bad weather in Arizona in December 2011. 
Figure[9^ shows the entire USNO light curve, illustrating 
that at least a few flux dips are present beyond the end 
of the CoRoT campaign. Figure [9 ]d overplots the IRAC 
and USNO data for the time window where we have some 
optical data that falls at or near the center of one of the 
strong IR flux dips. The ratio of the dip depth at 4.5 pm 
to the depth in the I band - A4.5/A/ - is about 0.3. In or- 
der to use this ratio to help constrain the occultation ge- 
ometry, we must make a correction for the fraction of the 
IR light that comes from the disk (assuming the disk flux 
is not significantly occulted). Based on a fit to the star’s 
SED, our estimate of f4.5(disk)/f4. 5 (star) is 1.3; hence 


the corrected dip depth ratio becomes A4.5/A7 ~ 0.7. 
Given the uncertainties, this is consistent with equal dip 
depths at both wavelengths, as one would expect from 
dust structures that are optically thick but cover only 
a small part of the star. If instead the dust structures 
are optically thin, this ratio would indicate an extreme 
deficit of small particles compared to the standard ISM 
mixture where the ratio should be ~0.10 (Indebetouw et 
al. 2005). 

4 . 4 . Mon-6975 

Mon-6975 is another essentially anonymous star that 
had not even been considered an NGC 2264 member un- 
til now. Our attention was drawn to this star by its 2011 
CoRoT light curve (it was not observed by CoRoT in 
2008). We have subsequently obtained optical and IR 
spectra of it from SOAR (see Appendix) which indicate 
a spectral type of M2. 5 and moderately strong Ha emis- 
sion (Ha equivalent width = 15 A). Good seeing images 
show it to be a nearly equal brightness optical binary 
with a separation of order 2.4" oriented NE-SW on the 
sky; the NE component has a nearly featureless spectrum 
that is much redder than that of the M2. 5 star. Despite 
the very red color, the NE component has Ha weakly in 
absorption, and no other spectral features that are de- 
tected with certainty. We therefore assume it is not a 
YSO and is only a line-of-sight companion to the SW 
component. Both objects contribute to the C 0 R 0 T light 
curve, though we assume that the flux dips arise only 
from the young star. We do not have a measure of its 
UV excess, nor a complete set of IRAC photometry in 
order to classify its disk. 

The C 0 R 0 T light curve for Mon-6975 shows a nearly 
constant “continuum” level and a large number of narrow 
flux dips with a r elatively wide range in both depth and 
width (Figure [lofc.). The periodicity of the flux dips for 
this star is not as obvious as for Mon-21 and Mon-1165, 
but an auto- correlation function yields a period of 2.8 
days. To guide the eye, we have placed a cross below the 
first strong dip in the fight curve, and then successive 
crosses spaced 2.8 days apart for the remainder of the 
light curve. Some of the flux dips for Mon-6975 are very 
narrow, with Gaussian er values ranging from 0.095 to 
0.13 days. Figure |TT7] < shows several of these dips, and 
the Gaussian fits wenave derived. 

We have an IRAC 3.6 pm light curve for Mon-6975 - 
it is overlaid on the CoRoT light curve in Figure [TT1 The 
IRAC light curve is similar in character to th eCJoRoT 
light curve in the sense that it shows a relatively flat con- 
tinuum level and two probable flux dips at the location of 
two of the strongest flux dips in the C 0 R 0 T light curve. 
The IRAC light curve is, however, much noisier, making 
an accurate assessment of the depth of the two flux dips 
at IRAC wavelengths difficult. However, it is clear that 
the derived ratio, after correction for the contribution to 
the IRAC flux from the disk, would be greater than 0.5 
and possibly consistent with 1.0. Therefore, as for Mon 
1165, the dip depth ratio is either consistent with op- 
tically thick occulting structures or with optically thin 
“clouds” whose dust mixture is extremely depleted in 
small grains. 

4 . 5 . Mon-1131 
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Fig. 8. — a) IRAC light curve for Mon-1165; 0.3 mag has been added to the 3.6 fi m (Ch. 1) photometry in order for the 3.6 and 4.5 
fim (Ch. 2) light curves to approximately align vertically. The plus signs and diamonds mark the two components of each dip pair; b) 
Cutouts from the IRAC light curve - with the 3.6 and 4.5 micron points shifted vertically to align with each other - at the times of the 
four well-sampled flux dip pairs; c and d) Illustration that the narrow flux dips for Mon- 1165 are reasonably Gaussian in shape. In this 
case, the y-axis for [4.5] has been both shifted and stretched, so that the dip has the same apparent amplitude in both channels. 


Mon-1131, also known as QT Mon, has CoRoT data 
from both 2008 and 2011. The 2008 light curve for this 
star is very different from the 2011 light curve, and is 
probably dominated by variable accretion - this light 
curve is shown and discussed in the Appendix. The 2011 
light curve (Figure |~i~2|: >) fits within the short-duration 
flux dip class, although compared to the other promi- 
nent members of the class, its flux dips are relatively 
broad. The dip period is 5.08 days. The two narrow- 
est dips have Gaussian a values of 0.27 and 0.30 days - 
their profiles are shown in Figure [12)3. The IRAC light 
curve for Mon-1131 is essentially uncorrelated with the 
CoRoT light curve. We also obtained a SOAR optical 
spectrum of Mon-1131. Based on measures of the depths 
of the TiO bands in that spectrum, we estimate its spec- 


tral type at Ml. 5. The Furusz et al. (2006) Multiple- 
Mirror Telescope (MMT) high resolution spectrum for 
Mon-1131 shows it to have a broad and structured Ha 
profile. Both its Ha equivalent width (>30A) and its 
UV-excess (Figure [23|) identify it as the most strongly 
accreting star amongst the short-duration flux dip group. 


4.6. Mon-1580 

Mon-1580 is another relatively little-studied member of 
NGC 2264. There is no published spectral type, however 
it is a known Ha emission line YSO (Reipurth et al. 1996; 
Furesz et al. 2006). Based on the Furesz et al. spectrum, 
we estimate the spectral type as Ml. We have a 2008 
CoR.o Flight curve, but no corresponding light cur ve f rom 
the 2011 campaign. The 2008 light curve (Figure [13^.) is 


16 


-0.4 ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' r 

. Mon- 1165 (2011) USNO 


-0.2 - . * I 



5900 5920 5940 5960 5980 

MJD-50000. (days) 


Mon— 1165 



o Ch2 


0.4 - b) 


5904 5906 5908 5910 5912 5914 

MJD - 50000. (days) 


Fig. 9. — a) USNO /-band light curve for Mon-1165; flux dips with depths up to about 0.4 mag appear to be present, b) Overlay of the 
IRAC and USNO light curve for Mon-1165, with relatively good data in all three bands for the flux dip near MJD 55907. 
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Fig. 10. — a) CoRoT light curve for Mon-6975. The plus sign under the flux dip at MJD 55906 marks the centroid of that dip. All of the 
other plus signs are simply placed at exactly 2.8 day intervals with respect to the first one; b) Gaussian fits to four of the narrow flux dips 
for Mon-6975 marked in the previous figure. All are well approximated by either a single Gaussian shape or the blend of two Gaussians. 
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Fig. 11. — Overlay of the CoRoT and IRAC light curves for Mon-6975; a zero-point shift has been added to the IRAC magnitudes so 
they approximately align with the CoRoT data. The two strongest CoRoT flux dips have counterparts in the IRAC light curve, and there 
is moderate correlation in the structure of the rest of the two light curves. 
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Fig. 12. — a) CoRoT light curve for Mon-1131; b) Illustration that some of the flux dips for Mon-1131 are reasonably Gaussian in shape. 


reminiscent of the 2008 Mon-21 light curve, in the sense 
that it appears to be the superposition of a spotted star 
light curve and a number of short duration flux dips. The 
spot light curve period and the flux dip period are the 
same ( P ~7.1 d ays) to within our ability to measure the 
periods. Figure [l3p shows the Gaussian fits to the two 
strongest flux dips. 

4.7. Mon-378 

Mon-378 is another system where the light curve is a 
superposition of a sp otted star light curve and narrow 
flux dips (Figure [l4p.). Mon-378 is a CTTS with an in- 
verse P Cygni Ha profile (Furesz et al. 2006), spectral 
type K5.5 (Dahm & Simon 2005), and location corre- 
sponding to a Class II YSO in the IRAC color-color di- 
agram. The CoRoT light curve is from 2011 - it was 
not observed in the 2008 CoRoT campaign. The IRAC 
light curve shows no flux dips and little correlation with 
the CoRo Plight curve. The Co Ro Plight curve shows a 
large amplitude spotted star waveform with a period of 
about 11 days and with a single prominent, narrow flux 
dip at MJD ~ 55919. However, after creating a model of 
the spotted star light curve as done for Mon-21, and sub- 
tracting it from the original light curve, a sec ond promi- 
nent, narrow flux dip becomes evident ( Figure [T4|3). This 
by itself would be a weak case that the dips are periodic. 
However, we also have a 2010 ground-based optical light 
curve, which shows (albeit poorly sampled) three appar- 
ently n arro w flux dips with the same ~11 day period 
(Figure [Tl}:). 

4.8. Mon-1076 

Mon-1076 (aka V338 Mon) is somewhat similar to 
Mon-378 and Mon-1580. It has a 2008 Co Ro Plight curve 


with low amplitude variations which may be due to star- 
spots but which are not stable enough to yield an obvi- 
ously periodic (sinusoidal) signature. The CoRoT light 
curve shows three prominent narrow flux dips at MJD 
54539, 54551 and 54556.5, consistent with a period of 
~5.8 days (Figure 15 %). Fits to the two best observed 
flux dips are shown m Figure The published spec- 
tral type for Mon-1076 is Ml(Dahm & Simon 2005). 
Mon-1076 was observed by CoRoT in 2011, but bad sys- 
tematics in the light curve preclude using it for any anal- 
ysis. 


4.9. Mon-566 

Mon-566 is a relatively late-type CTTS member of 
NGC 2264, with a published spectral type of M3. 5 
(Dahm & Simon 2005) and an Ha equivalent width of 
19.4 A(Dahm & Simon 2005). Its location in the IRAC 
color-color diagram places it in the Class II YSO cate- 
gory. It was observed by CoRoT in both 2008 and 2011. 
Mon-566 is quite faint in the optical, and therefore the 
CoRoT light curves are very noisy which makes inter- 
pretation of the light curves difficult. Nevertheless, both 
the 2008 and 2011 light curves appear to show periodic, 
narrow, shallow flux dips. The flux dip period is 4.5 days. 

Two other sets of data give us confidence that the 
CoRoT flux dips are real. First, a Spitzer 2011 light curve 
also shows low S/N flux dips of similar width and depth 
at the location of some of the CoRoT flux dips. Sec- 
ond, we have a USNO light curve for Mon-566 in 2011, 
and that light curve also shows flux dips coincident with 
those present in the CoRoT and Spitzer light curves. All 
three light curves are shown in Figure [16| 

The flux dips in the USNO light curve seem to be sig- 
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Fig. 13. — a) Co Ro T light curve for Mon-1580. This is another star where there appears to be a sine- wave-like, spotted-star light curve, 
superposed on which are narrow flux dips which recur for at least two or three epochs. Two such recurring dips are marked with diamonds 
and circles in the plot; b) Illustration that some of the narrow flux dips for Mon-1580 can be well fit by single Gaussians or blends of two 
Gaussians (dips A and C in Table |3|). Though not well-determined due to the relative shortness of the observing campaign, the period for 
the semi-sinusoidal (spotted-star) waveform and the average interval between dip systems (blue diamonds or red circles) is about 7.1 days. 
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Fig. 14. — a) Observed CoRoT 2011 light curve for Mon-378; b) Result of subtracting a model spotted star light curve, generated as in 
Sec. 3.1.1 for Mon-21, from the as-observed light curve; c) USNO light curve for Mon-378 from observations in 2010. The period for both 
the spotted-star waveform and the short-duration flux dips is about 11 days. 
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Fig. 15. — a) CoRoT light curve for Mon-1076. The first diamond is centered at the measured location of the first flux dip; the other 
diamonds are centered at the time predicted assuming a period of 5.80 days; b) The two most prominent flux dips and our fits for Mon-1076 
(dips A and B in Table [ 5 ]). 
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Fig. 16. — a) Overlay of the IRAC and CoRoT light curve for Mon-566. The magnitude zero points are shifted in order to align the 
two IRAC light curves and displace them above the CoRoT light curve for ease of comparison. The large black filled dots mark flux dips 
present in the CoRoT light curve; the magenta open circles mark flux dips that are present in the IRAC light curves, b) USNO 7-band 
light curve for Mon-566; flux dips with depths up to about 0.15 mag are present, coincident in time with some of the flux dips shown in 
the CoRoT/ Spitzer plot. The green diamonds are placed at the same times as for panel (a); the red diamond marks a dip beyond the end 
of the CoRoT campaign, where it would appear given the measured 4.5 day period of the other flux dips. 


nificantly deeper than those in the CoRoT light curves. 
We suspect this is because Mon-566 is very faint for 
CoRoT ( r = 17.7), and in some cases the background 
subtraction routine that is used in the CoRoT pipeline 
results in too small (or large) a correction. For Mon-566, 
the correction is probably too small and the apparent 
flux dip depths are therefore also too small. Most of the 
stars in our sample are considerably brighter in the opti- 
cal, so we do not believe their flux dips are significantly 
affected by imprecision in the background correction for 
their Co Ro Flight curves. 

5. YSOS SHARING SOME TRAITS WITH THE 
SHORT-DURATION FLUX DIP STARS 

Are there stars with light curves exhibiting superposed 
narrow and broad, periodic flux dips? Are there stars 
which have just one or a couple narrow- flux dips over 
the timescale of the CoRoT short runs? Do Weak-lined 
T Tauri stars (WTTS) ever show narrow flux dips in their 
CoRoT light curves? These questions could shed light 
on both the physical mechanism for the phenomena and 
the requirements on our line-of-sight to the circumstellar 
disks in these system. 

To answer the above questions, we critically examined 
the CoRo T light curves for all ~600 NGC 2264 members 
(CTTS, WTTS and unknown classification) in our cata- 
log for which we have CoRo T light curves as well as for 
600 field stars with CoRo T light curves obtained simul- 
taneously with the NGC 2264 members. The member 
stars included all of the YSOs identified in McGinnis et 
al. (2015) as AA Tau analogs. We find that two of the 


stars with AA Tau-type light curves do exhibit periodic, 
narrow flux dips in addition to the broad flux dips re- 
quired to be members of the AA Tau class. Figure [17] 
shows these light curves, with the red diamonds marking 
the location of the periodic, narrow flux dips. 

In each case, the period inferred from the narrow flux 
dips is equal (within the uncertainties) with that derived 
from the broad flux dips. Gaussian fits to the narrow flux 
dips yield widths that are in the same range as those mea- 
sured for the narrow flux dip stars of Table [T] For these 
composite light curves, the narrow flux dips always occur 
during the portion of the light curve near the minimum 
in light for the broad-flux dips. The amplitudes for the 
narrow flux dips vary considerably from epoch to epoch, 
but do not seem to correlate with the amplitudes of the 
broad flux dips. 

We have also identified another type of composite light 
curve - light curves with periodic features that are some- 
times narrow and sometimes broad (defined arbitrarily 
as FWHM/Period > 0.3). There are at least four such 
objects in our NGC 2264 dataset - Mon-441, Mon-928, 
Mon-379, and Mon-14132; all are included in McGinnis 
et al. (2015) as AA Tau analogs. Their light curves are 
shown in Figure [l8j It could be argued that these four 
stars should be added to our group of short-duration, pe- 
riodic flux dip stars. We prefer not to do that because 
we prefer to restrict the group to stars where all the dips 
in a CoRoT campaign are short duration. However, in- 
cluding these stars in the group would not change any of 
our conclusions. 
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Fig. 17. — View of a portion of the CoRoT light curve for two AA Tau analogs, Mon-250 and Mon-811, highlighting in each light curve 
several narrow flux dips that are spaced in time at the period defined by the broad-flux dips. The red diamonds mark the location of the 
narrow flux dips. The periods are P = 8.75 day and 8.0 day, respectively. 


Based on the fact that the periodic narrow and broad 
flux dips can be present in the same star and during 
the same campaign, the most significant deduction one 
can make is that the requirements for producing these 
two types of dips overlap. Because the AA Tau, warped 
disk model requires a specific viewing angle to the disk, 
that same viewing angle must also be compatible with 
production of the short-duration, periodic flux dips. 

Our visual re-examination of the other CoRoT light 
curves yielded two additional results. First, none of the 
stars identified as WTTS (nor any of the non-member 
stars) showed any narrow or broad flux dips whose mor- 
phology mimicked those of the AA Tau light curves or 
those of the periodic, short-duration flux dip stars. This 
is, of course, as one would expect if such dips are due 
to close-in circumstellar dust, which WTTS and non- 
members should lack. Second, we identified a set of nine 
stars all CTTS - whose light curves showed only a 
single, or in some cases a few, short-duration flux dips 
(but with no hint of periodicity). Light-curve snippets 
for a sample of these flux dips are shown in Figure [19} 
The widths and depths of these dips fall within the range 
found for the stars of Table [l] The physical properties 
(UV excesses; IR excesses; spectral types) of the stars 
with these single dips are not obviously different from 
the other types of variable extinction stars. We assume 
that the same physical mechanism(s) for producing the 
flux dips applies to this set of stars, but in their case 
other dust raising events either do not lift the dust high 
enough to occult our line of sight or yield dust opac- 
ities insufficient to yield a dip detectable in our data. 
However, when considering the statistics of what frac- 
tion of CTTS show any evidence of line-of-sight extinc- 


tion events, these objects should also be included in the 
calculation. A table of basic information for these stars 
(Table |5j) is provided in the Appendix. 

6. PHYSICAL PROPERTIES OF THE SHORT DURATION 
FLUX DIP CLASS 

Why do the CTTS in Table [T] show periodic, short- 
duration, often Gaussian shapea flux dips? If there are 
common physical characteristics that are shared by this 
set of nine YSOs, perhaps that would help allow us to as- 
sign a physical mechanism to the flux dips. We begin by 
examining the inferences we can make regarding the ori- 
entation of the circumstellar disks to our line of sight. We 
then discuss a variety of color-color and color-magnitude 
diagrams (CMDs) which sort the NGC 2264 YSOs by 
mass or accretion rate or IR excess. For the inferred 
properties in Table [l] and for the calculations in this sec- 
tion, we adopt a distance of 760 pc for NGC 2264 (Sung 
et al. 1997). The Sung et al. distance is in good agree- 
ment with a recent VLBI distance for a maser source 
in the cluster (d = 738+gQ pc; Kamezaki et al. 2014) 
and with a recent analysis of data for a newly-identified 
NGC 2264 pre-main sequence (PMS) eclipsing binary ( d 
= 756±96 pc; Gillen et al. 2014). 

6.1. Our View Angle to The Disks of the Narrow 
Flux-Dip Stars 

Practically the simplest, most straight-forward, test we 
can make to constrain the origin of the material causing 
the flux dips is to determine the inclination of our line of 
sight to the circumstellar disk for these stars. In order 
for a YSO to show flux dips from disk material passing 
through our line of sight, most models predict that the 
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Fig. 18. — Light curves for four NGC 2264 CTTS whose light curves show periodic flux dips that are sometimes broad and sometimes 
narrow. The red diamonds mark the location of the periodic, usually narrow flux dips - in each case, placing the first diamond near the 
bottom of the first dip, then spacing the other diamonds by multiples of the period we have determined. We, and McGinnis et al. (2015) 
include these stars with the classic AA Tau group, and they are shown as red circles with inscribed red plus signs in Figure [3] though their 
sometimes narrow flux dips suggests a link with the stars of Table [l] 


inclination of the disk to our line of sight must be * > 65° . 
To increase the sample size for this test, we include not 
only the stars in Table[l]but also the other stars discussed 
in §5 that have both narrow and broad flux dips. In total, 
we have eleven stars for which we have all of the relevant 
information (luminosity, stellar T e g, usinz, and rotation 
period). The input data, and our estimates of the line of 
sight angles are provided in Table 0} 

The periods listed in Tableware those derived from the 
flux dips. In several cases (Mon-21, Mon-56, Mon-378, 
Mon-1580), the light curves also show semi-sinusoidal 
variations plausibly due to photospheric spots, and in 
those cases the spot periods are indeed essentially the 
same as the periods derived from the dips (and thus 
consistent with “disk locking”). We assume for the 


other stars that their flux dip periods are also equal to 
their photospheric rotation periods. The v sin i values 
are either from VLT-FLAMES or Keck-HIRES spectra 
we have obtained for this program, or from reanalysis 
of MMT-Hectochelle spectra obtained by Furesz et al. 
(2006) . A description of the spectral analysis is provided 
in McGinnis et al. (2015). The stellar luminosities are de- 
rived from the measured J-band magnitudes, extinctions 
based on fits to multicolor photometry, the adopted dis- 
tance of 760 pc, and bolometric corrections from Pecaut 
& Mamajek (2013); the derivation of these luminosities is 
described in Venuti et al. (2014). Effective temperatures 
are derived from spectral types - either from the liter- 
ature or from our own spectra as described in the text 
and the appendix - and a conversion from spectral type 


24 




13.6 

j ^ 

13, 


13.8 

1577 | 1 

15.9 

5930 5935 

A* lit 


16 

r inf? 

16.1 

1 1 1 
CO 

o 

16.2 

— i — i — >— 


5910 5920 




15.6 

15.7 

15.8 


vA /ft 

: V 1 

15.2 

15.3 

m 

16 

16.05 

ih. - 

: 

i{ 

11559'j 1 

15.4 

11568 j]| ' 

16.1 

i - 

6991 : 

'1 ■■ ■■ 1 ■ ■■' 

15.5 


16.15 

■ 1 ■ ■ ■ ■ 1 ■ ■ - 


4540 4545 5920 5925 5930 5935 

MJD - 50000. (days) 


Fig. 19. — Isolated, short-duration flux dips found in the light curves of NGC 2264 CTTS, for stars with only one or a few such events 
per CoRoT campaign - and hence for stars not previously included in our variable extinction classes. Biographical data for these stars are 
provided in the Appendix (Table [5j . 


to T e ff appropriate for YSOs. In order to illustrate the 
impact of the temperature scale, we show results for two 
representative temperature scales in Table [4j one from 
Table 6 of Pecaut & Mamajek (2013; PM13) and one 
from Cohen & Kuhi (1979; CK79). 

For many of the stars, the calculation returns a value 
of sin* > 1.0, and in these cases we tabulate an incli- 
nation of 90°. A v isual summary of the results is pro- 
vided in Figure |20[ The mean values and uncertainties 
are: < sin* > = 1.16 ± 0.55 (CK79) and < sin* > = 
1.08 ± 0.53 (PM13) when including all eleven stars; and 
< sin i > = 1.03 ± 0.35 (CK79) and sin* = 0.96 ± 0.35 
(PM13) if one excludes the highly discrepant values for 
Mon-928. These estimates support the hypothesis that 
we are viewing these systems nearly edge on to their 
disks, and that the material passing through our line of 
sight to produce the flux dips plausibly lies within the 
disk or near the disk surface. It is clear, however, that 
individual sin* estimates have large uncertainties, mak- 
ing this result less robust than we had hoped. Assuming 
that T e ff and v sin * uncertainties dominate, and that they 
contribute approximately equally, the sin * estimates sug- 
gest the true uncertainties in T e g are of order 10% and 
in usin* are of order 20%. 

6.2. Stellar and Disk Properties of the Narrow Flux-Dip 

Stars 

Color-magnitude or color-color diagrams could provide 
clues to the physical mechanisms driving the photomet- 
ric variability in our three light curve classes (AA Tau 
analogs, stars with aperiodic flux dips, and the new 
short-duration flux dip class). If a mechanism depended 
on the stellar mass or luminosity, this should be reflected 
by the location of the class members in a CMD. If the 
mechanisms depended sensitively on our view angle to 


the circumstellar disk, then the mean Ay to the system 
as estimated from a J — H vs. H — K s diagram could 
differ amongst the three classes. If the mechanism(s) de- 
pend sensitively on the radius of the inner disk wall or 
the degree of disk flaring, that could be reflected in the 
IRAC [3.6]— [4.5] vs. [5.8] — [8.0] diagram. 

Figure [21~[ i compares the location in an optical CMD 
of the short duration flux dip stars to the set of classical 
AA Tau analogs (except as noted in footnote 22) and 
the aperiodic extinctor candidates in NGC 2264 as iden- 
tified in McGinnis et al. (2015). No reddening correction 
has been attempted. The primary conclusion we draw 
from this diagram is that the AA Tau analogs and stars 
with aperiodic flux dips appear to be bluer and brighter 
(higher mass), on average, than the narrow dip stars, 
with only a couple having g — i colors corresponding to 
M dwarfs. Secondarily, we see no obvious age difference 
between the three groups of stars or of these stars relative 
to the entire set of NGC 2264 members. As illustrated in 
Figure [22j the spectral types for the short-duration flux 
dip stars are also significantly later, on average, than 
those for the NGC 2264 stars we classify as “classic” AA 
Tau analogs, lending support to the belief that they are 
generally lower in mass. Based on a Monte Carlo simu- 
lation, there is less than a 2% chance the two groups are 
drawn fr om the same parent spectral type distribution. 

Figure [2T}> provides the J — H vs. H — K s diagram for 
the same set of stars, as a means to assess their red- 
dening. Since all of the variable extinction stars are 
CTTS, we assume that the Meyer et al. (1997) CTTS 
locus represents their intrinsic (unreddened) colors, at 
least for low mass CTTS. The aperiodic dip star located 
near H — K s = 0.4, J — H = 0.6 (Mon-325) is a high 
mass star (spectral type G6), and presumably dereddens 
back to the ascending portion of the main-sequence locus. 
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Fig. 20. — Histograms of the derived sinz estimates for the stars exhibiting periodic short-duration flux dips. The red dashed histogram 
corresponds to the estimates made using the CK79 temperature scale; the blue solid histogram corresponds to the calculations made with 
the PM13 young star temperature scale (their Table 6). 


TABLE 4 

View Angles to the Disk Axis 


Mon ID 

Period 

(days) 

v sinz 
(km/s) 

Luminosity 

(L©) 

Spectral 

Type 

Ref a 

Tett 

(PM13) 

T e g 

(CK79) 

Rstar(CK79) 

(R©) 

sinz(CK79) 

(R©) 

i-PM13 

(degrees) 

i-CK79 

(degrees) 

CSIMon-000021 

3.20 

19.2 ± 1.5 

0.957 

K5 

R02 

4140.0 

4395.0 

1.68 

0.719 

39.53 

45.97 

CSIMon-000056 

5.83 

11.1 ± 2.2 

0.774 

K5 

DS05 

4140.0 

4395.0 

1.52 

0.842 

48.20 

57.35 

CSIMon-000378 

11.1 

10.4 ± 0.4 

0.896 

K5.5 

DS05 

4080.0 

4295.0 

1.71 

1.333 

90.00 

90.00 

CSIMon-001076 

5.85 

9.8 ± 2.4 

0.256 

Ml 

DS05 

3630.0 

3681.0 

1.23 

0.908 

61.98 

65.20 

CSIMon-001131 

5.18 

8.8 ± 1.7 

0.545 

Ml. 5: 

S14 

3560.0 

3590.0 

1.91 

0.470 

27.60 

28.00 

CSIMon-001165 

5.50 

16.8 ± 2.3 

0.130 

M3: 

S14 

3360.0 

3360.0 

1.06 

1.711 

90.00 

90.00 

CSIMon-001580 

7.10 

11.4 ± 2.0 

0.265 

Ml: 

S14 

3630.0 

3681.0 

1.26 

1.260 

90.00 

90.00 

CSIMon-000250 

8.70 

9.8 ± 2.2 

1.229 

K3 

DS05 

4550.0 

4775.0 

1.62 

1.037 

70.29 

90.00 

CSIMon-000379 

3.69 

25.0 ± 1.2 

2.121 

K2 

DS05 

4760.0 

4955.0 

1.97 

0.920 

58.06 

66.86 

CSIMon-000811 

7.88 

13.5 ± 1.7 

1.076 

K6 

DS05 

4020.0 

4200.0 

1.96 

1.070 

78.50 

90.00 

CSIMon-000928 

9.92 

20.1 ± 2.8 

0.576 

M0 

DS05 

3770.0 

3915.0 

1.65 

2.490 

90.00 

90.00 


a Spectral type references: R02 - Rebull et al. (2002); DS05 = Dahm Sz Simon (2005); S14 - this paper 


Most of the variable extinction stars lie near the Meyer et 
al. CTTS locus, and therefore have little additional red- 
dening as inferred from their mean 2MASS photometry. 
One exception is Mon-6975, whose IR colors of H — K s 
= 0.42, J — H = 0.90 indicates relatively large inferred 
reddening. Mon-6975 is also the reddest variable extinc- 
tion star in the optical CMD (Figure 21 1 ). We conclude 
that outside of the flux dips, the other YSOs with nar- 
row and broad flux dips have at most modest reddening, 
and therefore that reddening corrections would not sig- 
nificantly alter their locations in optical CMDs such as 


Figure [Till . 

Figure |23^ shows the u — g vs. g — r color-color dia- 
gram (Venuti et al. 2014; Stauffer et al. 2014), discussed 
extensively in Venuti et al. (2014) as an accretion diag- 
nostic. YSOs lacking active accretion (WTTS) form a 
roughly inverted V-shaped locus in this plane, with late 
F stars located near g — r = 0.6, u — g = 1.6 and later 
spectral types becoming redder until about M0 at g — r 
— 1.4, u — g = 2.8. Later type M dwarfs maintain g — r 
~ 1.4 but have bluer u — g colors for later M subtypes, 
presumably due to active chromospheres. Accretion hot 
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Fig. 21. — Comparison of the location of members of the three different classes of variable extinction YSOs in an optical CMD (panel 
a) and in the near-IR J — H vs. H — K s diagram (panel b). The magenta curve in panel (b) is the 2MASS main-sequence locus (Pecaut 
& Mamajek 2013); the magenta dashed line is the Meyer et al. (1997) CTTS extension to the main sequence locus (transformed to the 
2MASS filter system). The small black dots are all probable NGC 2264 members (Cody et al. 2014), regardless of whether we have CoRoT 
light curves or not. A number of objects are marked both as broad dip and aperiodic dip members because their light curve in 2008 falls 
in one class and their 2011 light curve falls in the other class. 


spots primarily affect the u — g colors, displacing the star 
towards the bottom of the plot (bluer u — g colors) . The 
diagram shows that most of the variable extinction stars 
have weak UV excesses. The short duration flux dip stars 
have locations in this diagram indistinguishable from the 
members of the two other variable extinction light curve 
classes. 

Figure [23) o provides the IRAC color-color diagram, 
which allows one to separate stars with IR SEDs char- 
acteristic of Class I, II and III YSOs (Allen et al. 2004). 
The diagram shows that a significant fraction of the AA 
Tau analogs and aperiodic flux-dip stars have IR colors 
intermediate between Class II and III, whereas the short 
duration flux dip sample all have relatively red colors 
consistent with the Class II locus. The red IR colors 
imply that warm dust near the inner disk wall often pro- 
vides the majority of the flux in the IRAC bands (see 
column 7 of Table [l]). That, plus the expected lower se- 
lective extinction ax 4 /jm compared to the optical, may 
help explain why, when we do have simultaneous optical 
and IR light curves, we usually see only hints of the flux 
dips in the IRAC data. 

7. TOWARDS A PHYSICAL EXPLANATION FOR THE 
SHORT-DURATION FLUX DIPS 

In the previous sections, we have attempted to estab- 
lish quantitative data characterizing the short duration 
flux dips and the stars whose light curves exhibit these 
features. We summarize these properties: 


• About 5% (8 of 159) of the CTTS in NGC 2264 
show this type of light curve (periodic, short- 
duration flux dips) at least some of the time. Af- 
ter removing two objects we classify here as short- 
duration flux dip stars, McGinnis et al. (2015) con- 
sider 31 of 159 of the CTTS with CoRoT light 
curves as either AA Tau analogs or aperiodic ex- 
tinctors, corresponding to a 19% occurrence rate. 

• Another nine CTTS show one or a few short- 
duration flux dips in their CoRoT light curves. 
Therefore, the total fraction of CTTS with some 
evidence for flux dips due to dust structures inter- 
secting our line of sight is 48/159 or 30%. 

• For the stars of Table [I] the flux dip durations are 
typically a < 0.25 days, with a range from about 
0.10 to 0.4 days, based on Gaussian fits. 

• Many of the flux dips are well-fit by single Gaus- 
sian profiles, or in some cases, blends of 2 or 3 
Gaussians. 

• The depths are never large - the largest we have 
measured is of order 25%; most are only of order 
10 %. 

• The periods exhibited range from about 3 to about 
10 days. In some cases, we also can derive stellar 
rotation periods from superposed senri-sinusoidal 
flux variations in the light curves, and when we do 
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Fig. 22. — Spectral types for the stars in Table [T| compared to those for the classical AA Tau stars in NGC 2264 (same stars as in 
Figure pi. For one of the narrow-dip stars and six of the AA Tau analogs, there is no published spectral type and we have adopted the 
T e jj values from Venuti et al. (2014) to estimate spectral types. The spectral bins used were K1-K3; K3-K5; K5-K7; K7-M0; M0-M1.5; 
M1.5-M3; and M3-M4.5. 


so, the periods are equal or nearly equal to the dip 
periods. 

• There are significant changes in both depth and 
shape of the dips from epoch to epoch; in most 
cases the widths are more constant than the depths. 

• The flux dips persist for typically many orbital pe- 
riods. In at least two cases, periodic narrow dips 
are present in both 2008 and 2011. 

• In three cases (Mon-21 in 2008, first and last dips; 
Mon-56, several dips; and Mon-1165, all dips), 
there are paired flux dips with the dips separated 
well enough so as not to be “blended”. When this 
occurs, the separation in time between the two dip 
centroids stays approximately constant. 

• Where we have estimates of the inclination of the 
disk to our line of sight, that inclination indicates 
that we usually view the disk close to edge on. 

• Compared to the other CTTS in NGC 2264 with 
variable extinction light curves, those with short 
duration flux dips are, on average, later in spec- 
tral type, with the majority being early to mid M 
dwarfs. 

• In the IR.AC color-color diagram, the short- 
duration flux dip stars have colors which place 


them in the middle of the Class II distribution, 
whereas the classic AA Tau and aperiodic flux dip 
stars often have mid-IR colors in the transition re- 
gion between Class II and Class III. 

• A small number of stars show both broad (AA Tau- 
like) and narrow periodic flux dips; another small 
set show one or two narrow flux dips, with no evi- 
dence of periodicity. 

• For the very small number of events where we have 
sufficient coverage at both optical and IR wave- 
lengths, the ratio of the flux dip depth in the optical 
to that at 4 /im indicates that either the occulting 
structures are optically thick, or if the “cloud” is 
optically thin, then the dust grain mixture must 
be very depleted of small grains compared to the 
standard ISM composition. 

7.1. Inferred Characteristics Which Any Model Must 
Incorporate 

We assume that any physical model to explain the 
short duration flux dips will have as its starting point 
that the dips are due to our line of sight being inter- 
cepted by structures in or near the inner circumstellar 
disk that rotate through and occult our line of sight. A 
fundamental parameter of such structures is their physi- 
cal size. We can estimate the average size scale involved 
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Fig. 23. — Location of the short-duration flux dip class members in two diagrams sensitive to on-going accretion and warm dust. The 
dashed line in the left-hand plot separates stars that are strongly accreting (those below the line) from those with weak or no on-going 
accretion (Venuti et al. 2014). The box in the right-hand plot is the locus of Class II YSOs from Allen et al. (2004). Class III YSOs have 
colors near 0,0 in this diagram. The small black dots are all NGC 2264 members and candidate members, regardless of whether we have 
CoRoT light curves or not. 


from the parameters in Tableland Figure[3] From Fig- 
ure [3j the average value of the ratio of the dip FWF1M 
and the period for the short-duration flux dips is about 
0.15. Assuming the dips are due to material orbiting 
the star, that ratio should correspond approximately to 
the angular size of the structure as seen from the star, in 
this case indicating an angular size of about 0.95 radians. 
The average co-rotation radius for the stars in Tabic [I] is 
11.3 R.0. Assuming that the dust structures are located 
at or near the co-rotation radius, the average physical 
size of the structures would then be about 11 R 0 . For 
our stars in Table [I] the stellar radius R star ~ 1.5 R 0 , 
and hence the physical size of the structure corresponds 
to about 7 Rstar, or to about 0.05 AU. 

Further, the shallowness of the dips - usually just 5- 
10% - implies either that only a small part of the star is 
being occulted, or that the average optical depth through 
the dust structure is small. The lack of a flat bottom to 
the flux dips requires that the dust structure not be very 
small relative to the star, nor much bigger than the star 
(completely occulting it) and homogenous. 

Finally, we emphasize again the frequent occurrence 
of flux dips whose shapes are approximately Gaussian. 
In combination with the preceding two points, this char- 
acteristic places significant constraints on the geometric 
structure of the occulting object. 

With these considerations in mind, we offer four possi- 
ble physical models which could produce short duration, 
periodic flux dips in the light curves of CTTS. 


7.2. Model 1: The Short- Duration Periodic Flux Dips 
Fall within the Standard Warped Disk Model for 
AA Tau Light Curves 

Classical AA Tau light curves usually show consider- 
able structure within their broad flux dips, which is of- 
ten highly variable from epoch to epoch. These char- 
acteristics presumably attest to physical structures in or 
near the inner disk wall whose size, shape, and/or optical 
depth also change considerably on the timescales of the 
Keplerian orbital period at that radius. The duration of 
the AA Tau flux dips is set by the fraction of the orbital 
longitude range where the upper edge of the inner wall 
intersects our line of sight. Even if the height of the up- 
per lip of the disk wall varies smoothly with longitude, 
there will be certain lines of sight where the longitudinal 
range which intersects our line of sight is small, result- 
ing in a short-duration, periodic flux dip. More probably, 
the height of the inner disk wall will vary irregularly both 
in time and in longitude, with the highest points being 
potential sources for short-duration flux dips given an 
appropriate view angle to the disk. 

Arguments in favor of this interpretation of the data 
are that it is the simplest explanation, and it can fairly 
naturally explain stars with periodic flux dips having a 
range of (dip FWHM)/period. The simple mathematical 
formula proposed by Bouvier et al. (1999) results natu- 
rally in there being a width-depth correlation for flux 
dips produced by disk warps of varying inclination to 
our line of sight - as illustrated in Figure [24j Given the 
Bouvier et al. formula to describe the disk warp for AA 


29 


Tau, flux dips with maximum depths > 0.3 magnitudes 
and FWHM/Period > 0.25 (classic AA Taus) would be 
predicted to be somewhat more common than those with 
shallow but detectable depths, again as observed. The 
simple model also yields dip shapes that are reasonably 
close to Gaussian. 

In detail, however, some aspects of the data do not 
seem to fit with this simple model. The distribution 
of variable extinction stars in Figure [3] is somewhat bi- 
rnodal, with a deficit of objects with FWHM/Period ~ 
0.3 and Depth ~ 0.3, perhaps more suggestive of two 
populations rather than a single population. The ten- 
dency for the classic AA Taus to have earlier spectral 
types and to be bluer and brighter in the optical CMD 
argue against there being just a geometrical (view-angle) 
difference between the classic AA Taus and the short- 
duration flux dip sample. Instead, there could be a mass 
or temperature effect. For example, this could be ex- 
plained by differing magnetic field geometries as a func- 
tion of mass, over the K1 to M4 range in spectral types 
for our stars. Also, while the simple geometrical model 
can produce Gaussian flux dip shapes, the actual flux 
dip shapes for the classic AA Taus are highly structured, 
requiring the material in the disk warp to also have a 
highly non-uniform and variable distribution. If the disk 
warp model were adjusted to yield structured dips to 
match the classic A A Taus, it is not obvious it would 
yield predicted narrow flux dips that would often still be 
well-approximated as Gaussian in shape. 

7.3. Model 2: Flux Dips Produced by the Footpoint of A 
Funnel-Flow Accretion Column or by 
Inhomogeneities in a Dusty Disk Wind 

Romanova et al. (2004; 2008) predict that for YSOs 
with relatively low accretion rates, the accretion may oc- 
cur in stable, relatively long-lived “funnel flows”. Typi- 
cally, two tongues of gas arise from the inner disk, sepa- 
rated by about 180°. The gas originates from a relatively 
broad range in longitude and then gathers together into 
a narrower stream (it “funnels”). Some empirical evi- 
dence for these types of accretion curtains has been pro- 
vided through analysis of CTTS emission line profiles 
and their variability (Symington et al. 2005; Bouvier et 
al. 2006). If the gas flow originates from the region near 
the dust sublimation radius, and that radius is beyond 
the co-rotation radius, then it will start with a lower 
angular rotation rate than the star and a trailing spiral 
pattern will result, as illustrated in Figure 1 of Romanova 
& Kurosawa (2014). Dust could be entrained in this flow 
and survive briefly as the gas/dust mixture rises above 
the disk plane (though the dust will soon sublimate as 
it flows towards the star). Alternatively, Bans & Konigl 
(2012) have advocated that disk winds originating from 
beyond the dust sublimation radius may levitate suffi- 
cient dust to contribute a large fraction of the 3 to 4 
micron flux in YSOs, and if that wind were significantly 
non-axisymmetric it could be the source of potentially 
short-duration flux dips. 

We have constructed a toy model of an accretion cur- 
tain in order to illustrate constraints that the narrow dip 
profiles place on such a curtain. We tune the model to 
the flux dips in Mon-21. Specifically, the model must 
produce approximately Gaussian shaped flux dips with 
a ~ 0.2 days (FWHM ~ 0.45 days) for a system with 


-Rstar = 1-7 Rq, P = 3.2 days, and R co = 9.7 Rq. We 
place the curtain at the co-rotation radius. Based on 
the geometrical arguments from §7.1, we adopt a curtain 
width of 8 -R s tar- We divide the curtain into 8 vertical 
panels, each of base width 1 7? sta r (the left panel of Fig- 


ure 25 illustrates the geometry of the toy model). The 
panels are actually triangles, all meeting at a common 
apex at height //convergent above the midpoint of the cur- 
tain. For very large //convergent) the panel is essentially a 
vertical wall of width 8 Rst&r- For //convergent = 3 As tar - 
the curtain becomes an equilateral triangle. We adopt a 
uniform optical depth of 0.1 for dust within the curtain. 
Finally, we place the bottom of the curtain so that it is 
aligned with the bottom of the star as viewed from our 
vantage point. We place the edge of the curtain near the 
leading limb of the star, and let orbital motion (with the 
Mon-21 period) carry the curtain in front of the star. 

With the above idealized geometry, we wish to deter- 
mine if flux dips like that observed can be reproduced 
- and if so, what range of converge nt p oint heights is 
allowed? The right panel of Figure [25 shows the flux 
dips produced for values of the convergent height from 
20 -R s tar (bottom) to 1 i? s tar (upper blue curve). The 
model with a convergent height of 1 R s tar produces an 
approximately Gaussian shaped profile whose FWHM ~ 
0.4 days, within the range of dip widths observed for 
Mon-21. These results emphasize two things: (a) in or- 
der to produce a flux dip as broad as observed, an ac- 
cretion curtain must be quite wide (several times the 
stellar diameter) at its base; and (b) in order to produce 
the Gaussian shape, assuming a uniform dust opacity, 
the curtain must “funnel” down to a much smaller re- 
gion over a very small scale height. The alternative to 
this would be to assume instead that the curtain has 
an opacity distribution smoothly varying with longitude, 
peaking near its midpoint. 

What types of limits on opacity variations within the 
curtain could we detect? We constructed one model, for 
convergent height of 1 R s tar , where we set the opacity 
of one of the panels to half the normal value. The top 
model flux dip in Figure [25) > shows this result. For many 
of the flux dips illustrated in this paper, the observed dip 
shape is more symmetric than this model and we could 
exclude an opacity fluctuation of this order. However, 
some of the flux dips we have shown do have asymmetries 
which could perhaps be indicative of opacity fluctuations 
(rather than as a blend of two Gaussian waveforms). 

One criticism of this model is that it could simply be 
regarded as a reformulation of the first model. That is, 
the tip of a warped disk is the place where one might 
expect an accretion curtain to form, and whether one la- 
bels the dust structure there as part of the disk or the 
base of the accretion flow is, perhaps, only semantics. 
One reason to keep the two ideas separate is that if the 
line of sight inclinations to Mon-21 and especially Mon- 
1131 (Table [4]) are confirmed at a value near our cur- 
rent estimates (i ~ 40 and 27 degrees, respectively), this 
could not be accomodated within a warped disk model 
but could be compatible with the accretion curtain model 
if dust could survive in the accretion flow to heights of 
several stellar radii above the disk plane. In this case, the 
duration of the flux dip would be set by the inclination 
of the accretion flow relative to the orbital axis. 

The, on average, later spectral types and dustier disks 


30 


sc 

cd 

S 


n 

0 


1*>^v 

• * • • • 




•yr*r r 


••••• 


• • 


• •• • 

,• •••• 

••• 

• • • 


• •. . 


- 0.4 - 0.2 0 0.2 0.4 

Phase 


Fig. 24. — Model flux dips for a disk warp of fixed height and extent, where the inclination of our line of sight is varied from 65 (yellow) 
to 75 (light blue, at bottom) degrees in steps of 1.25 degrees. The model used is that from Bouvier et al. (1999), with the disk parameters 
that were appropriate to AA Tau - hmax = 2.7 i? s tar, inner disk radius Ri = 8.8 -Rstar, and 4> = d=180°. 


for the short-duration flux dip stars are not an obvious 
feature of the model but could possibly be made to fit. 
The lower luminosity of the star and the dustier disk 
could correspond to allowing the base of the accretion 
column to be located beyond the dust sublimation ra- 
dius, whereas the higher luminosity, less dusty disks of 
the classical AA Taus might have their accretion columns 
be dust free. Alternatively, the less dusty nature of the 
A A Tau group’s disks may simply reflect their higher 
stellar mass and the expected dependence of disk evo- 
lution timescales on stellar mass (longer timescales for 
lower mass stars). 

7.4. Model 3: The Short-Duration Periodic Flux Dips 
Originate from Dust Associated with Vortex 
Instabilities Located Near the Inner Disk Rim 

The short durations and small amplitudes of the flux 
dips could arise from mountains on the surface of the disk 
- that is, features standing up above their surroundings. 
Such features would either be low-relief, sweeping across 
only the bottom of the star’s face but optically thick, or 
taller, but diffuse and of moderate optical depth. Moun- 
tains arising from the disk would naturally cross the star 
once per local orbit, while varying from one crossing to 
the next due to the internal turbulence. The issue is how 
to lift the mountains above the plains. 

One way to produce vertical relief on the disk surface is 
through vortices driven by Rossby wave or baroclinic in- 
stabilities. Sub-millimeter and millimeter interferomet- 
ric imaging of nearby transition-disk YSOs has revealed 


several stars with large azimuthal asymmetries in surface 
brightness at separations from the star of tens of AU (van 
der Marel et al. 2013; Isella et al. 2013; Perez et al. 2014). 
The asymmetries are interpreted as millimeter-sized par- 
ticles concentrated in long-lived hydrodynamical vortices 
rotating in the disk plane (Lyra & Mac Low 2012; Lyra 
& Lin 2013). Similar vortices occurring within a fraction 
of an AU of the star might account for the periodic flux 
dips we observe. The vortices have higher surface densi- 
ties than their surroundings and thus could increase the 
column along lines of sight to the star in some range of 
viewing angles. However to account for all the narrow 
dips, the vortices would need to obscure the star in 5- 
10% of all systems, and thus to increase the height of the 
starlight unit-optical-depth surface near the dust subli- 
mation front by 5-10% of the radius - or more if not 
every disk has a vortex near the relevant radius. 

In the existing hydrodynamical models, the surface 
density enhancements at the vortices are about a fac- 
tor of two (Lyra et al. 2009; Meheut et al. 2010; Lyra 
& Mac Low 2012; Meheut et al. 2012). We construct 
a model consisting of a star of 0.5 M 0 and 2 Rq, with 
effective temperature 4000 K. Orbiting the star is a disk 
accreting at 1O _8 M 0 yr -1 . The disk’s surface density 
is 1000 g cm^ 2 at 0.1 AU and falls off inversely with 
the radius. An inner edge is created by rolling off the 
surface density inside 0.1 AU as a Gaussian with half- 
width 0.015 AU. The disk’s outer edge is defined by an 
exponential cutoff with radial scale 40 AU. To prevent 
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Fig. 25.— Illustration of the geometry used by the toy model (left panel), and light curves produced by allowing orbital motion to pass 
the model accretion curtain through our line of sight to the star (right panel). The curves plotted are for convergent heights of 20 R a tar, 2 
fistar, 1.5 -Ratar, and 1 R s tar- The top (magenta) curve is also for a convergent height of 1 Rstar, but with one panel of the curtain set to 
half the opacity of the other panels. 


the outer disk front obscuring the star before the inner 
edge as the system is tilted over to higher inclinations, we 
give the outer disk a lower dust abundance. The dust- 
to-gas mass ratio is 1% inside 0.2 AU and falls off as 
a Gaussian with half- width 0.2 AU toward dust-to-gas 
mass ratio 10 -4 . The disk material’s opacity comes from 
well-mixed compact spherical grains chosen to match the 
observed spectra of protostellar sources (Preibisch et al. 
1993). We place a vortex at 0.1 AU, giving it a surface 
density contrast of two and profiles in the radial and az- 
imuthal directions that are Gaussians with full-widths at 
half maximum of one-quarter the radius and 60° , respec- 
tively. 

We translate the disk’s surface densities into volume 
densities using a temperature profile taken from the mid- 
plane of a Monte Carlo radiative transfer calculation in a 
similar but axially-symmetric density distribution, lack- 
ing the vortex. This places the model disk in approxi- 
mate hydrostatic equilibrium. We then carry out a 3-D 
transfer calculation to obtain the distributions of tem- 
perature and radiation intensity in the disk with vor- 
tex, using the Bjorkman & Wood (2001) algorithm and 
Lucy’s (1999) method as implemented by Turner et al. 
(2012). Isotropic scattering is included. We solve the 
transfer equation on a grid of rays (Yorke 1986) includ- 
ing the scattering and thermal emission source terms, 
to construct synthetic images of the system as viewed 
at different points around the vortex’s orbit. We con- 
duct photometry on the images and assemble them into 
a lightcurve. The resulting V-band photometric ampli- 
tude is 1% or less at all viewing angles. The model vortex 
thus has too little vertical relief to explain the dips’ am- 


plitudes. 

To account for the possibility that vortices can even- 
tually reach contrasts up to an order of magnitude rela- 
tive to their surroundings, as suggested by long-duration 
2-D hydrodynamical calculations (Regaly et al. 2012), 
we carry out a further run in which the surface density 
contrast is increased to ten. The V-band amplitude ex- 
ceeds 3% when the system is inclined between 81 and 
84°, peaking at 7% for inclination 83°. Therefore, even 
this seemingly extreme vortex cannot easily explain ei- 
ther the full range of the narrow dip amplitudes or the 
number of objects in which they occur. However, vor- 
tices concentrate pebble- to boulder-sized rubble (Barge 
& Sommeria 1995), and collisions between the rubble 
particles could yield fragments that raise the abundance 
of sub-micron dust even more than by a factor of 10. De- 
termining whether the increase makes vortices a viable 
explanation for the narrow flux dips will require better 
models of the dust sources, sinks and transport. 

7.5. Model 4 : Flux Dips Produced by A Spiral Arm in 
the Inner Disk 

A planet embedded in a disk exerts a gravitational 
force that raises a spiral arm in the surrounding ma- 
terial. The pattern co-rotates with the planet and is 
strongest near the planet’s location (Rafikov 2002). In 
particular, a planet orbiting interior to the disk’s inner 
rim could raise a bump on the rim that rotates with the 
planet’s orbital period, which is shorter than the period 
at the rim. Hot Jupiters orbiting mature stars cluster 
near periods 3-5 days and eccentricity zero, while less 
massive planets show no such clustering. Hot gas giants 
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cannot easily grow in situ, as reviewed by Helled et al. 
(2013). Instead they probably formed outside the snow 
line and reached their present orbits either by migrat- 
ing through the disks that gave them birth, or by later 
dynamical interactions with third bodies (Baruteau et 
al. 2013). Detecting planets interacting with the disks 
would help separate these possibilities. 

To test whether planet-disk interactions could cause 
the narrow dips, we construct a model using the same 
0.5 Mq star and disk as in the vortex picture above. The 
vortex is replaced by a 5 M e planet orbiting at 0.07 AU, 
which raises a spiral overdensity in the disk following the 
shape of Muto et al. (2012) eq. 1, which is derived from a 
linear analysis by Rafikov (2002) and therefore is strictly 
valid only for planets with masses less than about 1 . 

We build the density and temperature distribution using 
the same procedure as for the vortex model, and again 
conduct synthetic observations using the radiative trans- 
fer code. The results are shown in Figure [26| where the 
head of the spiral arm passing in front of the star pro- 
duces a narrow dip with optical amplitude 2.3%, in the 
range observed among the stars in NGC 2264. A more 
massive disk or higher dust abundance would likely yield 
larger amplitudes. Also, the spiral arm’s star-facing side 
is hotter than the axisymmetric model used to set up the 
approximate hydrostatic equilibrium, so that restoring 
the balance would make the atmosphere expand, increas- 
ing the dip amplitude. However whether there is time to 
reach hydrostatic equilibrium is unclear, since the spiral 
arm propagates around the disk at the planet’s orbital 
frequency. A more serious difficulty with this picture 
is that it yields dips in less than 5% of all randomly- 
oriented systems, even if every star has a similar hot 
planet. The planet-driven bump simply does not extend 
high enough above the surrounding disk to obscure the 
star over a wide range of viewing angles. A more massive 
planet could yield a taller bump, but accurate modeling 
would require non-linear numerical hydrodynamical cal- 
culations. If the amplitude can be made large enough, 
then among the four pictures we have considered, the 
planet-raised spiral arms could best explain any individ- 
ual system in which the dips are regularly-spaced and 
repeat more often than the orbital period at the dust sub- 
limation radius. However, unless hot Jupiters are much 
more frequent in NGC 2264 than amongst field stars, this 
mechanism could not be invoked for the majority of the 
periodic, narrow, flux dip stars. 

8. CONCLUSIONS 

This paper, and the others in this series (Alencar et 
al. 2010; Cody et al. 2014; Stauffer et al2014; McGinnis 
et al. 2015) demonstrate how high cadence, long dura- 
tion, sensitive, photometrically stable light curves can be 
used to infer physical properties of YSOs that are oth- 
erwise impossible to investigate. In this paper, we have 
identified a new YSO light curve class characterized by 
short-duration, shallow, periodic flux dips. Depending 
on which physical mechanism is responsible for produc- 
ing these short-duration flux dips, detailed study of this 
set of stars has the potential to significantly improve our 
understanding of the structure of the inner disks and ac- 
cretion flows in YSOs and the role of instabilities and 
migrating planets in shaping the inner disks and their 
dust content. 


Regardless of the precise physical mechanism for pro- 
ducing the short-duration flux dips, if we are correct in 
assigning the semi-sinusoidal variations to photospheric 
spots and the short-duration flux dips to dust in or near 
the co-rotation radius, then our observations for at least 
Mon-21 (2008), Mon-378, and Mon-1580 confirm that 
their stellar rotation rate is locked to the Keplerian ro- 
tation period at the inner disk. This has been a long- 
standing prediction of some models for the angular mo- 
mentum evolution of YSOs (e.g., Konigl 1991; Collier- 
Cameron and Campbell 1993) but with generally only 
indirect evidence to support it. 

We offer here one possible interpretation of our data, 
which attempts to incorporate as many as possible of 
the constraints on the physical mechanisms we have de- 
scribed in Sec. 7. Assume that extinction events can be 
produced both by a warped inner disk and by dust en- 
trained within an accretion flow of the sort envisioned 
in the Romanova et al. models. Assume further that for 
the youngest, least evolved disks, our line of sight to the 
inner disk may be obscured by optically thick dust in the 
outer, flared disk. Adopt the idea proposed in Keane et 
al. (2014) that transition disks, as well as having par- 
tially cleared inner disks, also have less flared (or less 
dusty) outer disks. Given those thoughts, the stars with 
isolated, narrow flux dips (Figure 19) could be systems 
where dust in an accretion flow has survived to an unusu- 
ally high distance above the disk plane. We see only one 
or two of these dips per 40 day light curve because only 
rarely does the dust survive long enough to reach a height 
to intersect our line of sight. These systems are, there- 
fore, viewed at somewhat smaller inclinations than the 
other variable extinction stars. The periodic, narrow flux 
dip stars are systems viewed at larger inclinations where 
it is more likely the dust will generally be present in the 
accretion flows. Finally, the classic AA Tau analogs are 
the systems with more evolved less dusty - outer disks 
where we see into the inner disk and where a warped in- 
ner disk can then occult our line of sight. The systems 
with both broad and narrow flux dips have extinction 
events produced by dust in both the warped disk and in 
the accretion funnels. We offer this scenario primarily as 
a way to focus future studies. The main results of this 
paper are the empirical constraints on the narrow flux- 
dip light curves and their host stars provided in §6 and 
§7, which any model of these events must explain. 

A general concern with any model for the variable ex- 
tinction stars is simply that there seems to be too many 
of them in NGC 2264. There are roughly equal num- 
bers of stars with each of the three light curve types 
(AA Tau, aperiodic dip, and short-duration dip), with 
about 39 such stars in all, corresponding to about 25% 
of the total number of CTTS for which we have CoRoT 
lig ht cu rves. If one includes the single-dip stars from Fig- 
ure 19 then about 30% of the CTTS have some evidence 
for oust clumps within 0.1 AU of the star intersecting 
our line of sight. That seems high, given an expectation 
that the inner disk wall height should be of order 15% of 
the disk radius at that location (Dullemond et al. 2001; 
Muzerolle et al. 2003), even allowing for some YSOs hav- 
ing warped disks. A provocative solution to this concern 
is suggested by the claim from Klagyivik et al. (2013) 
that there are many more eclipsing binary members of 
NGC 2264 than should occur by chance if their axes were 


33 



4535 4540 4545 4550 4555 

MJD - 50 000 days 


Fig. 26. — Photometric signature of the spiral density wave launched in a T Tauri disk by a planet orbiting at 0.07 AU near the disk’s 
inner edge, according to a simple near-hydrostatic model. The synthetic images at top have red, green and blue channels corresponding to 
the V, K and 3.6 ft m bands. The image at left is 0.4 AU across and the close-up at right is magnified ten times. Both show the spiral 
wave’s peak about to cross the star’s face. The system is 82° from face-on. At bottom is the V-band lightcurve (blue curve) overlaid on the 
starspot-subtracted CoRoT lightcurve of Mon-21 (gray points). The synthetic lightcurve’s period has been adjusted to match the observed 
3.2-day dip interval, but the amplitude has been left unsealed. 


randomly oriented, from which they conclude that the 
EB axes are preferentially aligned and that we are near 
the plane of their orbital motions. If this is also true 
for the angular momentum vectors of the single stars in 
NGC 2264, then we would also detect a higher fraction of 
variable extinction stars than would be true for a random 
orientation of axes. Other authors have concluded that 
the rotational velocities (v sin i and periods) of stars in 
open clusters are compatible with a random orientation 
of axes and incompatible with any strong alignment of 
axes (e.g., Kraft 1970; Jackson & Jeffries 2013), so this 
either argues against such an effect in NGC 2264 or that 
NGC 2264 is unusual in this respect. 
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Facility: Spitzer (IRAC), CoRoT, CFHT (MegaCam), 
VLT (FLAMES). 


APPENDIX 

SEPARATING THE SPOTTED-STAR AND FLUX 

DIP SIGNATURES IN THE 2008 LIGHT CURVE 
FOR MON-21 

In §4.1.1, we used the Mon-21 light curve itself to gen- 
erate a model for the portion of the light curve we at- 
tribute to spots, and subtracted that model from the data 
in order to derive an estimate of the light curve struc- 
ture due to variable extinction. In order to demonstrate 
that this procedure is adequate to the task, we offer here 
another means to estimate the flux dip light curve. 

All modern studies agree that the dominant physical 
mechanism which drives the visual appearance of opti- 
cal light curves for most WTTS is the rotational modu- 
lation of star-spots that are distributed asymmetrically 
over the stellar surface. When these spots are at mod- 
erately high-latitude, which is typical, the light curve 
shapes can be approximately sinusoidal. Dozens of the 
NGC 2264 CoRoT light curves fall into this category. 
Given that large a library, it seemed reasonable to us 
that one of those light curves might have a shape that 
closely resembled the smoothly varying portion of the 
Mon-21 light curve. We believe that is indeed the case. 
Figure [27] shows the 2011 CoRoT light curve of Mon-103, 
a YSO with spectral type very similar to that for Mon-21 
(K6 vs. K5 for Mon-21). The period for Mon-103 is just 
very slightly greater than for Mon-21 - 3.3 days vs. 3.15 
days. The red curve in Figure [27] is the template wave- 
form for the Mon-21 model, simply shifted in zeropoint 
to match phase with Mon-103 and multiplied in ampli- 
tude to match also; the waveform shapes are indeed very 
similar. 

We can use the Mon-103 light curve as a template to re- 
move the spotted-star component from the Mon-21 light 
curve. In order to provide as good a fit as possible, we 
make the following slight modifications to the observed 
light curve: (a) we do a zero point shift in time to align 
the two beginning points in phase; (b) we shift the Mon- 
103 light curve counts by a multiplicative factor and a 
zero point offset so that the mean magnitude and aver- 
age amplitude match that of Mon-21; and (c) we multiply 
the x-axis timescale by 0.935 in order to correct for the 
slight period difference between the two stars. Figure [28[i 
overlays the two light curves, and Figure [28) ) shows the 
Mon-21 flux dip light curves that results from subtract- 
ing the modified Mon-103 light curve from the observed 
data. The main dips appear essentially the same as in 
our original procedure (Figure [IJ:); the dip widths and 
depths (Table[3]) differ by less than 10% between the two 
renditions. 


MON-1131 IN 2008 

In 2011, the CoRoT light curve for Mon-1131 shows 
flux dips that recur periodically at P ~ 5.1 day intervals. 
The light curve shows little other structure except the 
flux dips. The 2 008 light curve shape is much more com- 
plex (see Figure 29) and probably owes its shape to sev- 
eral physical mechanisms whose effects are superposed. 
Yet, the mean magnitude is nearly the same during the 
two epochs, and the full range of the variations is also 


nearly the same (about 0.3 mag). The 2008 light curve 
simply has very different character. In terms of the quan- 
titative classification scheme of Cody et al. (2014), the 
2008 light curve would be classified as stochastic (with 
M = 0.12 and Q = 0.71), meaning the variations are 
fairly symmetric about the mean and essentially aperi- 
odic. There is still some power in the auto-correlation 
function at P = 5.1 days, but other variability domi- 
nates. 

NEW SPECTRA AND SPECTRAL TYPES 

We obtained new spectra for several of the stars dis- 
cussed in this paper, primarily in order to provide spec- 
tral types or v sin i for stars for which those data were 
lacking in the published literature. In particular, we ob- 
tained three Keck HIRES spectra (all on the same night) 
for Mon-21, a single Keck HIRES spectrum for Mon- 
1165, SOAR Goodman spectrograph medium resolution 
spectra for Mon-1131 and Mon-6975, and an INT low- 
resolution spectrum of Mon- 1165.. 

Snippets from these spectra are shown in Figure [30] 
For Mon-21 and Mon- 1165, we used the Spectroscopy 
Made Easy (SME; Valenti & Piskunov 1996) and BIN- 
MAGS (Oleg Kochukhov; http://www.astro.uu.se/ 
oleg) packages to derive v sin i estimates, as described 
more completely in the companion paper by McGinnis 
et al. (2015). 

For the lower resolution, INT and SOAR spectra of 
Mon-1131, Mon-1165 and Mon-6975, visual inspection 
showed that all three stars were early to mid M dwarfs. 
We therefore used spectral indices based on the depths 
of the observed TiO bands to infer spectral types; more 
specifically, we used the relations in Stauffer (1982) and 
Stauffer & Hartmann (1986) for field M dwarfs and re- 
lations in Herczeg & Hillenbrand (2014) calibrated for 
YSOs. The derived spectral types are Ml. 5 for Mon- 
1131, M3 for Mon-1165, and M2. 5 for Mon-6975. For 
Mon-1131, the field star relation yields M2 and the YSO 
relation yields Ml - we take the average; for Mon-6975, 
both relations yield M2. 5. We consider these spectral 
types preliminary because the relations are based on flux 
calibrated spectra, and we did not have the data needed 
to flux calibrate the SOAR or INT spectra. 

COORDINATES AND COROT ID’S FOR STARS 
WITH SINGLE OR JUST A FEW 
SHORT-DURATION FLUX DIPS 

As discussed in §5, a small number of the YSOs in 
NGC 2264 for which we have CoRoT light curves have 
either just one or a couple short-duration, shallow flux 
dips in their light curves. They therefore do not fall 
into any of our previously defined variable extinction 
classes. These stars are listed here in Table 0 All of 
them are Class II systems based on their IRAC photom- 
etry, though in some cases their IR excesses are relatively 
small. 
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Fig. 27. — CoRoT 2011 light curve for Mon-103 ( CoRoT SRaOl 223980447), a YSO whose light curve shows a fairly stable, nearly 
sinusoidal pattern likely attributable to a relatively high latitude spot-group. The red curve is the template waveform derived from the 
Mon-21 2008 CoRoT light curve, as discussed in §4.1.1. 


TABLE 5 

Basic Information for YSOs with Isolated Flux Dips 


Mon ID a 

2MASS ID 

CoRoT* 

CoRoT* 

SpT d 

Ha EW (A) 

CSIMon-000119 

06412100+0933361 

223985987 

223985987 

K6 

10.6 

CSIMon-000342 

06405573+0946456 


616872592 

M4 

21.1 

CSIMon-000577 

06414382+0940500 


616895846 

K1 

5.9 

CSIMon-000650 

06410098+0932444 

223980807 

223980807 

K1 

6.4 

CSIMon-001038 

06402262+0949462 


602095739 

M0 

2.7 

CSIMon-001114 

06393339+0952017 

223957142 

223957142 

Ml. 5 

19.7 

CSIMon-001559 

06393276+0904137 

223956963 




CSIMon-001568 

06410684+0902161 

223982299 

223982299 



CSIMon-006991 

06392200+1006233 

223953770 

223953770 




a The Mon IDs are our internal naming scheme for stars in the field of NGC 2264 — see Cody et al. (2014). 
k CoRoT identification number in 2008 SRaOl campaign. 
c CoRoT identification number in 2011 SRa05 campaign. 

^ See Cody et al. (2014) for the sources of the spectral type and Ha equivalent width (EW) data. All of these are in emission. 
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Fig. 28. — a) The as-observed CoRoT light curve for Mon-21 overlaid by the 2011 light curve for Mon-103, after adjusting the latter light 
curve by zero point shifts and linear stretches in both axes; b) Light curve resulting from subtraction of the modified Mon- 103 light curve 
from the Mon-21 light curve. 
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Fig. 29. — CoRoT light curve for Mon-1131 in 2008. 
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Fig. 30. — Top panels show low-resolution spectra obtained with the Goodman Spectrograph at SOAR for two of the short-duration 
flux-dip stars. Both are moderately early M dwarfs (Ml. 5 for Mon-1131 and M2. 5 for Mon-6975). Bottom panels show Ha emission 
profiles present in Keck HIRES spectra for Mon-21 and Mon-1165. The Ha profiles are typical of modestly active T Tauri stars; the central 
absorption feature for Mon- 1165 may be partially or wholly attributable to imperfect sky subtraction. 
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